CHEMICAL ENGINEERING SCIENCE 
GENIE CHIMIQUE 


1959 








Mass transfer to a turbulent fluid with and without chemical reaction 


T. K. SHerwoop and J. M. Ryan 


Chemical Engineering Department, Massachusetts Institute of Technology, Cambridge, Mass. 
(Received 12 December 1958) 


Abstract 
mass transfer with a simultaneous, fast, irreversible reaction. Data on dissolution of a rotating 


\ new theory based on a turbulent boundary layer model is developed to describe 


cylinder of benzzoic acid into aqueous caustic are compared with this theory and with the 
theories based on the film and penetration models. 

Comparison of the data on dissolution into water with a form of the Deissler analogy points up 
the fact that the various modern analogies between heat, mass and momentum transfer predict 
that the Stanton number should be proportional to y/fat high Prandtl or Schmidt numbers, whereas 
the data show direct proportionality. This suggests a basic error in the usual procedure by which 
analogies are developed. 


Résumé 
pour deécrire le 


Les auteurs développent une théorie nouvelle basée sur la couche limite turbulente 
transfert de masse avec une réaction instantanée rapide et irréversible. Les 
résultats sur la dissolution d’un eyvlindre d’acide benzoique animé d’un mouvement de rotation dans 
une solution de soude caustique sont comparés avec cette théorie et avec les théories de films 
et de pénétration. 


La comparaison des résultats sur la dissolution dans eau au moyen des formules d’analogie 
de Deissler indique le fait que les analogies diverses entre le transfert de chaleur, de masse et de 
quantité de mouvement prévoient un nombre de Stanton proportionnel a 4/f pour de grands 
nombres de Prandtl et de Schmidt, tandis que les résultats montrent une proportionnalité directe. 
Ceci suggére une erreur de base dans le procédé usuel de développement des analogies. 


Zusammenfassung—Auf der Grundlage eines Modells der turbulenten Grenzschicht wird eine 
neue Theorie der Stoffiibertragung mit umkehrbarer Reaktion 
Die Werte der Auflésung eines rotierenden Zylinders aus Benzoesiure in wissrigen 


gleichzeitiver schneller, nicht 
abgeleitet. 
Laugen werden mit dieser Theorie und mit der Theorie der Film- und der Durchdringungs- 
vorstellung verglichen. 


Ein Vergleich der Auflésungswerte in Wasser mit der Analogie nach Deissler zeigt, dass 


die verschiedenen modernen Analogien zwischen Wirme-, Stoff und Impulsiibertragung eine 
Proportionalitét zwischen der Stanton—Zahl und y/f bei hohen Prandtl- oder Schmidt-Zahlen 
voraussagen, wihrend die Messwerte eine direkte Proportionalitaét zeigen. 


grundsiitzlichen Irrtum bei der tiblichen Ableitung der Analogien hin. 


Das deutet auf einen 


UNDERSTANDING of the process of mass transfer and chemical reaction, most of the theoretical 


between phases is basic to the development of 
such important chemical engineering operations 
as gas absorption, solvent extraction, distillation, 
drying and leaching. Although industrial opera- 


tions frequently involve simultaneous diffusion 


and experimental investigations of mass transfer 
between two phases have been concerned with 
the physies of the transfer of a single diffusing 
substance. 


In comparing theory with data on heat or mass 
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transfer to a turbulent stream, it has usually been 
found convenient to treat the case of a round tube 
through which the fluid passes in steady, fully- 
dev eloped turbulent flow. This 
about which much is known, but is complicated 
by the fact that the bulk stream temperatures or 
concentrations vary with distance from the tube 


is a situation 


inlet, so that analysis of data to obtain transfer 
coetlicients is often difficult. Use of a rotating 
evlinder avoids this difliculty, since the tempera- 
radial 


ture or concentration depend only on 


distance. 


EXPERIMENTAI 


Benzoic acid was melted, cast and machined 
to fill a 0-80 mm-deep recess, 63-5 mm-long, cut 


in a 30-5 mm-o.d. stainless or lacquered-steel 


eviinder. The benzoic acid was machined smooth 


and flush with the steel. The cylinder so prepared 
was rotated on a vertical axis at the centre of a 
10-2 em-i.d. glass container holding 1555 cm® of 
liquid. The recess holding the acid was at the 
vertical centre the 20.3 em-high 
Ball bearings were provided at both top and 
bottom, and the cylinder rotated by a variable- 
speed drive operable over a speed range from 
2 to 6000 rev, min. 
neoprene seals at top and bottom were held in 


of contaimer. 


Plexiglass end plates with 


place by external steel tie rods, and the whole 
was immersed in a constant-temperature water 
bath. 

Special pains were taken to exclude air, both 
to avoid the introduction of CO,, and to eliminate 
the air bubbles which tended to collect at the 
rotating surface. the container was 
through a 6mm hole in the top plate, sealed by 
a rubber membrane during a test. No free air 
liquid surface existed in this hole. A 0-40 mm-wide 
annulus, 25 mm long between the hole in the top 
plate and the rotating cylinder, formed a liquid 


Access to 


seal. 

A first series of tests was made with distilled 
CO,-free water, taken from a vessel kept in the 
constant-temperature bath. The tests varied in 
length from 0.5 to 4 hr, during which times the 
water became 10 — 40 per cent saturated in benzoic 
acid, and the diameter of the test cylinder, 
measured before and after the test, decreased by 
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about 2 per cent. Eight to ten samples of 5 — 6 em* 
each were taken at intervals during the test, using 
a hypodermic needle to withdraw them through 
a rubber membrane. Samples were immediately 
replaced by the same volumes of water. Twenty- 
five tests at 25°C, two at 33°C and one at 17°C 
were made over a range of Reynolds numbers 
from 136 to 108,400. 

of tests with NaOH 
solutions was carried out in similar fashion, but 
the tests were shorter in duration (14 min mini- 


A second set aqueous 


mum). Samples were neutralized with identical 
additions of standard HC] and back-titrated with 
the alkaline solution used in the test. Twenty-five 
tests were made at 25°C, using initial NaOH 
normalities from 0.0090 to 0-64; in two cases the 
initial solution contained both sodium hydroxide 
and sodium benzoate. 

DATA 


REDUCTION of} 


In the case of the dissolution tests using water, 
graphs of In(l — Cy, C4.) versus time gave 
k, was calculated from 
the measured slopes. Small corrections were 
made in the measured values of C 4, to allow for 
the replacement of the small samples by water. 

Similar corrections for sampling were made 
in the tests with caustic. In addition, allowance 
was made for the fact that &, varied during 
the test as NaOH concentration decreased. 
If k, varies linearly with the bulk NaOH con- 
centration Cp, the slope being k,; y/C4,, then 
1/yIn{1 — y [(Cyo)snitiar ~ Cro )/C ae } Should be 
linear in time. This was found to be the case, and 
k, at zero time (k,;) was calculated from the 
slopes of such graphs. Since the changes in NaOH 
concentration during a test were small, the result 
is quite insensitive to the assumed value of y, 


excellent straight lines ; 


and k,; is accurately established. 

Values of k, and k, obtained by these pro- 
cedures were calculated using an area based on 
the measured length and arithmetic mean of 
initial and final diameters of the cast test section. 
Surface speed u, was calculated from the mean 
diameter and the measured speed of rotation. 
The known physical properties of the solutions 
were employed in the calculation of Nx. 

SEIDELL’s [1] data on solubility of benzoic acid 
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in water and sodium benzoate were used; the 
water data were checked closely by measurements 
of the solubility of pieces of the cast test sections. 
Cuanc’s [2] values of the molecular diffusion 
coeflicient of benzoic acid in water were used. 
These are checked closely at 25°C by Kine and 
Bropie [8], but are 10 per cent higher than 
reported by Brennerr Lewis [4]. The 
Schmidt number for acid—water at 
25°C was taken as 870. 


a benzoic 


I 
Mass Transfer without Chemical Reaction 
(1) Results 


The data for twenty-eight tests in which benzoic 


PAR’ 


physical dissolution 


acid was dissolved in water are plotted on Fig. 1 
as k versus Ny. Twenty-five were at 25°C ; 


three at other temperatures are corrected to 25°C 


Uy 


by assuming k to vary with o °® 4, o being the 
Schmidt number. 


close to a well-defined line (to avoid confusion on 


The points fall reasonably 


the graph, this line is not shown). 
The most extensive comparable data are those 
of Ersenserc et al. [5], who carried out similar 


tests, employing several sizes of rotating cylinders 


and containers. The cylinder radius, R,, was 
varied from 0-97 to 2-99em and the container 
radius R,, from 3.03 to 6-84cem. Though the 


ratio of gap width to inner cylinder diameter 


varied from 0-104 to 4-88, an excellent correlation 
was obtained using Ny, based on R,,; the effect 
of R, and gap width appeared to be unimportant, 
at least for the geometries tested. o was varied 
from 835 to 11,500 by using different diffusing 
species, and Nx, found to be proportional to 
ao °*4, The solid line shown on Fig. 1 is a close 
their hundred or more 
870). It 
is evident that the agreement between the present 
data and those of E1isenpere et al. is remarkably 


representation of one 


tests with benzoic acid and water (¢ 


good. 
Also shown on Fig. 1 are five tests reported 


Kine and Bropre [3], who used a 2-2 


by 2-2 cm 
diameter cylinder of benzoic acid rotated in 
500 cm’ of water. The recent results of BENNETT 
and Lewis [4], also shown, have been adjusted 
to be proportional to 
1-0 cm, 


rotated in a 3 cm container. The four sets of data 


to 25°C by assuming & 


oO O44, 


Their cylinder diameter was 
would appear to fix the relationship for benzoic 
acid—water within 10 per cent to a high degree of 
probability. 


2) Theory physical dissolution 

As described in a recent review by one of the 
authors [6], the various analogies relating mass, 
heat and momentum transfer agree well with the 
considerable amount of data on mass and heat 


transfer for fully developed turbulent flow in 
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Data on rate of solution of cylinders of solid benzoic acid rotated in water. 
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round tubes. These are developed from the 


following rate equation : 
dC 


N (D €) (1) 
dy 


where N is the molal flux, C the concentration, 
y the distance from the phase boundary, and D 
diffusion 


and « are the molecular and eddy 


coefficients. D may be assumed constant in 
dilute solutions, but « varies from zero at the wall 
to very large values (relative to D)) in the main 
turbulent stream. 

In order to obtain integrated forms of equation 
(1), it is necessary to relate « and y. In many of 
the theories ¢ is assumed to be equal or proportional 
to the eddy viscosity, and so obtained from known 
velocity gradients. von KARMAN [7] employed 


the generalized uu’ y* correlation in smooth 
tubes for this purpose and obtained good agree- 
ment with data on heat and mass [8] transfer 


His 


integration of equation (1) assumed « to be zero 


in tubes for values of o from 0.5 to 25, 
for 0 < y* < 5. Kays and Brorkiunp [9] show 
a similar analysis to agree well with their data 
on heat transfer from a rotating cylinder to air. 
The theory fails completely, however, at large 
values of o. 

In liquid systems D is so small that the slight 
the 
important, and in this range the eddy viscosity 


turbulence very near wall can be quite 
cannot be obtained from existing data on velocity 
profiles. Various authors, therefore. hav ce proposed 
empirical relations relating « to y, and so obtained 
relations for the wall flux N.. The most successful 
analyses of this kind appear to be those of 
Deissier [10] and of Lin et al. [11]. 

For the 


writes : 


region of 0 < y* < 26, Detss_er 


ev nur y . exp | n? ut y*)] (2) 


where the empirical constant n is determined 
from data at high Schmidt numbers to be 0.124, 
and u* is related to y* by the generalized velocity 
correlation for smooth For uy <4, wt 
may be taken equal to y*, and for u* <1 the 
equation reduces to « vy = n‘(y*)*. Equation (2) 
was developed and tested for flow in round tubes, 
but will be assumed to apply very near the surface 


tubes. 


of a rotating cylinder, where most of the resistance 
to mass transfer at high Schmidt numbers is 
concentrated. In the region near the wall Lix 
et al. assume ¢« v to be proportional to (y*)*; 
the result is not greatly different if their form is 
used in place of Detssier’s. 


On 
equation (1) becomes : 


with N. assumed constant, 


u 


integration, 


d (y*, @) (3) 


Here € 
the wall, @ is introduced to represent y/(+,,/p). 


is the concentration at distance y* from 


and o is the dimensionless Schmidt group, v D. 
The function é(y*, ¢) is the indicated integral. 


mass transfer coetlicients 
(2) will be 
26, and the 
26 to the bulk 
The latter may be 


In order to obtain 
Stanton 
employed from the wall to y 


and numbers, equation 
Reynolds analogy used from y 
stream concentration at m,. 


written : 


/ dC du (4) 


whence : 


o(Y% ) — 6(26) = uy V(2,f) — 12-85 (5) 


Nog 


This increment varies from about 3 to 8 as Nx, 
10.000 to 100,000. Since it 
small part of the total integral, it will be taken 


increases from is a 
as constant at 5-75, equivalent to using equation 
(3) with the limit y* 200 in place of yy . 

From equation (3) it follows that the total 
integral may be written ; 


2 as 
keal2> Ng ale ¢ (200, o) (6) 


where 7,2, is replaced by (1 2)fpu2. The top 
line of Fig. 2 represents this integral and provides 
a basis for calculating k° ; it checks the DetssLer 
asymptotic equation for high values of co. 

It seems not unreasonable to assume that the 
conditions very near the wall should be the same 
for a rotating cylinder as for flow in a round tube 
at the same f. Away from the wall, however, the 
turbulence pattern is quite different in the two 
systems. 
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Fig. Graph of é(y*. @) 


Let R, be the radius of the rotating cylinder, 
and R, be the radius of the outer fixed cylinder 
confining the fluid in an annulus of gap width 

R,. Let the Reynolds number Nx, be 


c a 


“f} » » , , ‘re 
defined as 2K, u,v, where u, 


speed of the rotating cylinder, At 


is the peripheral 
very low 
speeds the flow is laminar, as in a narrow-gap 
rotating-cup viscometer. There is an intermediate 
speed range in which a series of three-dimensional 
vortices may be noted [12, 13, 14], above which 
the motion is turbulent. The present study is 
confined to the fully turbulent regime. 

Since the turning torque is the same at any 
radius, 

(7) 

Townsenxp [15] gives the following expression 
for the turbulent flow condition : 


d(u/r) 


v Ss 
dr pP r Uy ( ) 


TBe er 

, is the time-mean product of the 

fluctuating components of velocity in the radial 

and tangential directions. The velocity profile 

over the principal portion of the gap has been 

shown by Taytor [16] and Wenpr [17] to be 
given closely by : 


where v, u; 


ur = 0-5u, R, (9) 


If the eddy viscosity «, is defined by : 


d (u/r) 


dr (10) 


— T£.,,=(v+e)r 
and the friction coefficient is based on cylinder 
radius and cylinder peripheral velocity, then 
equations (9) and (10) may be combined to give : 


(11) 


This suggests that «, and « are constant except 
near the inner and outer walls of the rotating- 
cylinder system. 

Kquating «, and e, equations (2) and (11) may 
be employed to relate « to y* for the entire region 
from wall to bulk fluid [18]. This is not worth 
while for the analysis of data at high values of o, 
however, since the intersection of equations (2 
and (11) occurs at y* 47 for Np, = 750, 
y’ ~ 88 for Ny, — 11,500 11-3 for 
Na. 53,200, The resistance beyond y* 47 


is very small for o 


and y* 


870 (benzoic acid in water), 
so the use of equation (11) would hardly change 
the results of the analysis based on DeEIssLEer 
to y" 26 and the Reynolds analogy to y, . 

(3) Comparison of theory with data on physical 


dissolution 


In order to compare equation (6) with the data 
on physical dissolution, it is necessary to have 


values of f for rotating cylinders. These are 


obtainable from the work of THropoRsEN and 


Recier [19], who employed both smooth and 


rough cylinders of several sizes rotating in oils, 


kerosene, water and air. Their data are well 


represented by the single line A-A on Fig. 3. 








2Rwuw 


f = 
‘Qe. 


Fig. 3. Friction coeflicients for rotating cylinders. 
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That far s with gap width is suggested by the 
data of TayLor [16] and Wenpr [17], but the 
nature of the variation has not been established. 

Experimental and calculated values of Ng, are 
compared on Fig. 1, which shows lines representing 
equation (6) for both the Drissier and the Linx 
et al. functions. The comparison is seen to be poor, 
especially at low values of Ny,. where the dis- 
crepancy is about 50 per cent, 

Both theoretical lines are based on analyses of 
transfer processes in tubes, and it may be that a 
different 
employed for rotating cylinders. It is also possible 


relation between « and y* should be 
that Turoporsen and Reerer’s values of f are 
inapplicable, since their cylinders were rotated in 
a large space and all of the mass transfer data 


were obtained with confined cylinders. Fig. 3 


suggests that finite gap widths give steeper / 


curves, more in line with the mass transfer data 
of Fig. 1. 

A more interesting speculation, however, relates 
to the possibility that the form of equation (6) Is 


wrong. Following von KArMAN, « has been 


taken to be a function of y*. Suppose, instead, 
that the proper distance parameter were assumed 
to be the dimensionless group y**, detined by : 


Y Tw 


w ec YU I 


” 
— 


(12) 


U,gvp , 


i.e., ev = f, (y**), instead of ev — /, (y*). 


Dresser and Lin et al. compare their “ analogies” 
with heat and mass transfer data over quite 
of Np, (10,000 to 100,000). 
Assuming Dertss.er’s function /,(y*) to have 
50,000, then the use 
of f,(y**) would require the values of k°/u, as 


narrow ranges 


been well chosen for Na 


calculated by Derissier to be multiplied by the 
ratio of y/(f, 2) at Ny, to vy (f 2) at Ng, = 50,000. 
The latter value is given by Turoporsen and 
Recier as 4/0-00825. Adjustment 
procedure of the dashed line of Fig. 1 representing 
Lin et al. would result in a curve of predicted 
values checking the best line through the data 


by this 


points within 10 per cent. Similar adjustment of 
the Dertssier line would give similarly good 
agreement between predicted and experimental 


values. 


There are several reasons to believe that « v 
should be a function of y** instead of y*. Most 
important is the success of the many correlations 
of experimental mass and heat transfer data in 
the form of graphs of j,, or jp versus Ny,. Where 
the drag is due to skin friction, these follow 
f 2, not y/(f 2), as called for by equation (6). 

If « is taken to be a function of y*, as was done 
by von KArmMAN and later workers, the mass 
and heat transfer analogies take the form : 


Na = 77 J (ir 


is assumed to be a universal function of 


(13) 


If u* 
y** instead of y*, the result is similar except 
that y(2 f) becomes 2 / in the second term on 
the right. The Reynolds analogy still holds for 
o 1, 
variation of Ng, with o is unchanged, since the 


and at high values of o the nature of the 


first term on the right is then negligible. 

The velocity distribution 
u* and y* seems well established by the experi- 
ments of Nikurapse, Reicuarpr others. 
Most of the older data were obtained at high 


universal relating 


and 
Reynolds numbers, with little variation in (ff. 


Recent 
higher values of u* at any specified y*; use of 


data of Corcoran and Sacer [20] give 


y** instead of y* would introduce a correction 


factor in the right direction. 
ScHLICTING’s [21] review of the history of the 


u* ~y* correlation indicates that the choice of 
y* asa parameter was based largely on dimensional 
It would be of interest to explore 
in the velocity 


arguments, 
the possible advantage of y** 
correlation. 

Whether or not the modified distance parameter 
y** has application to velocity correlations, it 
may be concluded that « v for mass and heat 
transfer must be expressed as a function of wall 
distance in terms of y** rather than y* if Ng, is 
to be proportional to f at high values of e. 

Fig. 4 demonstrates the general nature of the 
Jv ~Ju ~Jf 2 correlation for rotating cylinders, 
and supports the statement that Ng, is pro- 
portional to f rather than y/f. The data of Fig. 1 
are represented by the line for the Eisen Bere et al. 
correlation at ¢ = 870. Several results of Correr 


and Scumiprt [22] for vaporization of naphthalene 
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from a 2-44-cm diameter cylinder rotating in air 
are seen to agree well, though o was less than 
1 per cent of the values for the liquid systems. 
The two sets of data [9, 23] on heat transfer to 
air fall somewhat low at high Reynolds numbers, 
but do not agree well with each other. The recent 
data of SeBan and Jounson [24] on heat transfer 
from a 6-35-cm cylinder rotating in oil and in 
water (for Ny, > 2000, with free 
convection) fall some 10 
mass-transfer data. The friction data of Turo- 
DORSEN and Recier are in reasonable agreement 
with the experimental curves for j, and jy. 
Fig. 4 represents a remarkable correlation, 
in that it brings together data on heat transfer, 
mass transfer and friction for cylinders rotating 
in gases and liquids over a range of Reynolds and 
Schmidt numbers of more than one-thousand 


fold. 


negligible 


15 per cent above the 


Part II 

Mass Transfer with Chemical Reaction 
(1) Theory 
The theoretical treatments of mass transfer 
with chemical reaction have generally led to 
expressions defining a multiplying factor to apply 
to the coeflicient for transfer without chemical 
reaction. These have been based on two models : 
(a) the stagnant film, and (b) an interface being 

continually replaced by fresh surface. 


rapid, irreversible reaction 


A 
ren? 
4 ‘ \ 


Comparison of data on heat transfer, mass transfer and friction for rotating cylinders 


The first, developed originally by Harta [25] 
on the basis of the stagnant-—film model, leads 
to the relation : 


keg Dp Cry 


ky , Dal te (14) 


This is for the case of a solute 4 entering the 
fluid and reacting irreversibly and instantaneously 
with solute B diffusing toward the phase boundary 
from the bulk fluid (see Fig. 5). The transfer 


Q 











—.— 


Fic. 5. Film model emploved by Harta. 

coefficient & , is defined as the ratio of the molal 
flux N , to the concentration difference (C 4. — 0) 
and k, is the corresponding transfer coefficient 
Dx, for the diffusion of A from interface to 
bulk without chemical reaction (no B present). 
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Retaining the stagnant film model, Surrwoop 
and We! [26] have shown how the Harra theory 
may be modified to allow for the dissociation of 
electrolytes in water, where the molecular diffusion 
coeflicients D, and Dy, have little meaning. For 
solution of benzoic acid in water, with dissociation 
of benzoic acid neglected, their result may be 
written (see Appendix) : 
ks 2RT x’ 
ke D, F? 


ls 


where F, is the Faraday and A_, A_ are the 
limiting (zero concentration) ionic conductances 
of OH~ and Bz 
concentration of 


Cro (» 1~* Ce " 
Li 1)| . J 0 (15) 


CABo ( ‘Aw 


ions, respectively. C4, is the 


acid at 


to 
benzoic the phase 
boundary ; Cy. and C4,, are the bulk solution 
concentrations of sodium hydroxide and total 
sodium, respectively. It follows from equation 
(15) that because of the rapidity of OH~ diffusion, 
k , increases some 50 per cent faster with Cp, C4. 
than suggested by the Harra form. 

The second model involves transient penetration 
of diffusing species into a renewable interface. 
The use of this model for the development of a 
theory of 
second-order reaction 
DancKwerts [27] and others [28, 29, 30, 31], 
without allowance for varying diffusion rates of 


transfer with simultaneous, rapid 


has been described by 


the different ions formed by dissociation, The 


calculated increase of k, with Cp. C4, is not 


X 
greatly different than predicted by Harra [32]. 

It is now proposed to employed a third model 

the turbulent boundary layer*. Following the 
pattern of the Harta theory, assume that 4 
reacts rapidly and irreversibly with 6 moles of B 
at some plane a distance y;, from the surface. 
Then : 


(C0 — °) (16) 


, q@ 
N, @ = (UR, o4) E. 


*Note added in proof: Three recent articles treat diffus- 
ion and chemical reaction in laminar boundars layers. 
These are: Frrepitanper S. K. and Lirr M., Chem. 
Engng. Sci. 1958 7 229; Porrer O. D., Trans. Inst. Chem. 
Engrs. 1958 36 415; CuHamere P. L. and Youne J. D., 
Phys. Fluids 1958 1 48. 


$ (Yo. op) —$(¥R. ep) (17) 


With no B present 
and : 


A diffuses without reacting, 


d (Y% ’ a4) (18) 
This assumes that the limit Cy corresponds to the 
limit ye, an assumption shown by Rannie [83] 
to be acceptable for ¢ < 1. 
ky  $(Yo +74) 
k, ¢ (Yr- 74) 
( ‘Bo d (Yo 
bc tee [d(yo ’ op) 


a4) 


(19) 
$(Y¥R+op)| 


In order to evaluate this relation, it is necessary 
to relate « to y* and obtain values of the integral 
represented by ¢. Before doing this, however, it 
is of interest to inspect three limiting cases, for 
which « need not be known. 


(a) oy =~ og. Equation (19) reduces to: 


hg C po l 
kg OC 4a {1 — [6 (yR)/o (ys )]} 
Cu 1 
= (20) 
bC ay [1 — (hak 4)] 
whence : 
A i Cro 


ky ” oc (21) 


fu 


which is identical with the Harra equation, 


since D, = Dy. 


(6) wR <1. ©,» approaches zero at the wall, 
and in the limit 4 (yg, ¢) 


equation (6). Then : 


>Up o, as seen from 


ka $ (Yo » 74) 
YR A 
, d (Yo ) a4) 
d (Yo ’ Cp) 
$(Yo +?) — YR oR 
YR A 


whence : 


88 
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h ky fo Cc 

A 8 By > (25) 
ky Ka log Wy 

If the coefficient k° is represented as being pro- 
portional to a power function of D(k, « D*), 


ta = (22) “fa 


A 


then: 


Da Chao 
D, OC Aw 


(26) 


The condition yz - 
k4/k, is greater than 20 at o, = 1; 
than 7 at o, = 10; greater than 3.5 at o, 
or greater than 1-7 at a, = 1,000. 


greater 
100, 


(c) kya/k,y near 1-0, 
being very large, so that the value of o is not 


wa.2. t. " — . ; 
rhis is the case of y, 


important, and the result is the same as for 
o4 Op: 

ha 14 C Bo 

Ks OC 44 
These three limiting situations detine the variation 
of k, within quite narrow limits, though no 
assumption as to the variation of the eddy 
diffusivity has been made. 
the DerissLer 
follows by substitution into 
equation (19) of the integral defined by equation 
(3). Computed values of this integral as a function 


The complete theory, 
expression for «, 


using 


of y* and o are given in Fig. 2; y, is taken to 
be 200 for the reasons given earlier. The ordinates 
of Fig. 2 represent the mass transfer resistance 
from the wall to any y*; the use of these values 
is seen to lead to a theory similar to Harra’s 
but with resistance per unit length varying with 
distance from the wall. 


) Results 


Fig. 6 shows the data obtained for the dissolu- 
tion of benzoic acid into aqueous NaOH solutions 
of several normalities. The range of variables 
(Nae Cgo) Was limited by the attainable precision 
and reproducibility. At low NaOH concentrations 
the effect of the caustic was small and the ratio 
k4/k4 sensitive; at the highest normalities and 
Reynolds numbers the tests were necessarily of 
short duration and suflicient benzoic acid was 
dissolved to distort the rotating cylinder. The 
points plotted, however, represent initial rates. 


with chemical reaction 


1 applies approximately if 


3. 6. Data on solution of benzoic acid into aqueous 
NaOH compared with several theories. 


The data follow the general trend of the several 
theories, and show the curvature predicted by 
the model. It 
however, that the points fall some 0 


turbulent transport is evident, 
50 per cent 
above the line representing equation (19), the 
discrepancy being the greatest near the middle of 
the The Harra equation fits the data 


somewhat better than the new theory ; since the 


range. 


film model is not realistic, this must be considered 
fortuitous. 

The lack of agreement between the data and the 
boundary layer theory is not attributable to the 
possible error in the analogy approach discussed 
in an earlier section, since here we are concerned 
only with the ratio of coefficients with and 
without chemical reaction. More probably, the 
difficulty lies in the failure of the theory to allow 
for the increased speed of diffusion of OH™~ ions 
in a multiple-ion system. The data points fall 
roughly half-way between the boundary layer 
theory based on molecular diffusion, and the film 
theory based on ion diffusion. They fall about as 
much above the theory for molecular diffusion in 
a boundary layer as the film ion-diffusion theory 
falls above the film molecular-diffusion theory. 
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Allowance for the increase in C ,,. with accumula- 
tion of sodium benzoate near the wall would affect 
both ordinate and abscissa of the points but 
would not greatly change their positions relative 
to the theoretical lines. 

It appears diflicult to extend the theory further 
until more is known about turbulent mixing. In 
the zone where the reactants meet, there is a 
fine-grained turbulence of small eddies. Some of 
these are acid and some basic, though the region 
may be statistically neutral. The lines of mean 
concentration, as represented by Fig. 5, probably 
zero concentration. 


cross rather than meet at 


Gross mixing is accomplished by the motion of 


the eddies; ultimate mixing necessary for chemical 
The 


development of a complete theory must await the 


reaction by molecular and ion diffusion. 


solution of the problem basic to most studies 


involving turbulence — the mechanism of the 


dissipation of eddies. 


APPENDIX 

Sueawoop and Wet (26) emploved average concentration 
gradients and concentrations in using the Vrwocrap and 
McBatn equations for ion diffusion. It is now noted that 
this is unnecessary ; integrated forms are not difficult to 
obtain in certain cases. 

For example, the second case treated by Suerwoop and 
Wer (non-dissociating acid diffusing and reacting with 
base) involves the diffusion of Bz~, OH and Na* 
in the right-hand portion of the film pictured in Fig. 5 


hots 


The Na®* flux is set equal to zero 
RT A, 
Pr 


A, G, 
G, c, 
| A.C, 


ALG 

- (28) 
AC 
where the subscripts refer to anion and cation, and the 
prime to OH-. G represents the gradient dC dy. The 
condition of electrical neutrality requires the C, be equal 
to the sum of C_ and C This is introduced and the 
equation integrated to obtain : 


, | c-'| 
o/h vo a 


The constant is evaluated in terms of the concentrations 
at the two boundaries. The equation for N’, similar to 
equation (28), is written and C_ eliminated to give : 


2A. dC, 
=a. Fae i | > 


(29) 


constant 


RT 2. 


“53 (31) 
a 


where : D 


Integration and substitution of limits leads to equation 
(15). This is identical with the result obtained by Suer- 
woop and Wri 
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NOTATION 


a — exponent in equation (26) 

b =moles of B reacting with one mole of A 

C = concentration g mole, cm 

C4q Cag = concentration of reacting solutes A or 
B in bulk fluid stream yg mole 
concentration of diffusing solute A at 


gz mole 


en? 


em 
3 


phase boundary 
concentration in bulk stream 
concentration of cation, anion 
diffusion coeflicients of 
em’, sec 
96,500 Cg equiv. 
dimensionless 


dC dy 


g mole em 
g equiv, em* 
molecular 
species A, RB 
Faraday, 
friction factor 
concentration gradient 
mass-force conversion factor 

unity in ¢.g.s. system 
Ns Ns.2* (for mass transfer) 
Nst Np,” 3 (for heat transfer) 
mass transfer coeflicient, Ne 
mass transfer coeflicient of species A 
or B for transfer without chemical 


C Aa cm sec 


reaction 
initial (zero time) mass transfer coefli- 
cient in batch test em sec 
mass transfer flux g mole (sec) (em?) 
empirical constant in equation (2 
Prandtl! number 
Reynolds number 
Schmidt number 
Stanton number 
mass transfer flux at wall 
ion flux g equiv (see) (em*) 
V Twhe P = Ug VU 2), OF ty y (f 2) 
radial distance from axis em 
gas constant 8-314 J (°K) (g mole) 
= radius of container em 
radius of cylinder em 
temperature *K 
time em sec 
dimensionless velocity u/Q 
value of u* corresponding to bulk 
stream velocity 
= average velocity of main stream, (rela- 
tive to wall) cm /sec 


2Ry yp * 
» D 


K thy 


g mole, (sec) (em?) 


mean velocity 





[32)} 
[33]} 


Mass transfer to a turbulent fluid with and without chemical reaction 


fluctuating velocity, tangential em? sec = eddy conductivity cm? /sec 
peripheral speed of cylinder surface — cm/sec , A’ = limiting (zero concentration) ionic 
fluctuating velocity, radial cm /sec conductances of anions (A) (em? /)(V) (g equiv.) 
film thickness em limiting ionic conductance of cation 
distance from wall or phase boundary em (A) (em*) /(V) (g equiv.) 
dimensionless distance from wall yQv viscosity g /(sec) (em) 
dimensionless distance defined by density g/cm 
equation (12) Schmidt number, Ng; or Prandtl 
value of y* corresponding to bulk number, Np, 
stream concentration shear stress at radius r dyn /cm?* 
> value of y* at reaction zone shear stress at wall dyn cm? 
rate of change of ky with Cp, during » «= kinematic viscosity, = p/p cm? /sec 
batch test ¢(y*, o) = integral defined by equaton (3) 
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Zusammenfassung 
Teilchen von deren Form untersucht. 
Stromungsbereich (Re 


ks wird die Abhangigkeit der Sinkgeschwindigkeit isometrischer (regulirer) 
Fiir den Stokesschen (0 < Re - 
- 1000) werden mit Hilfe des von Wape.t. eingefiihrten Kugelwertes 


0, 2) und den Newtonschen 


der isometrischen Teilchen formelmissige Zusammenhinge zwischen der Sinkgeschwindigkeit 


und der Teilchenform aufgestellt 


Abstract 
their shape is 
(0 < Re 0.2) and Newtonian (Re 


investigated. \ 


mathematical 


The relationship between the settling velocity of isometric (regular) particles and 


equation, embracing the regions of Stoke 


1000) types of flow, between the above qualities is put 


forward in terms of sphericity as defined by Wapr.t 


Résumé 
triques (réguliéres) et leur forme 
et de Newton (Re 


Dans les domaines d'écoulement de Stokes (0 < Re .- 


Les auteurs ont etudice la relation entre la vitesse de sedimentation des particules isome- 


- 1000) on a établi sous forme de formules les relations entre la vitesse de 


sedimentation et la forme des particules a Taide de la constante introduite pour des sphéres par 


WaADeLt. 


1 EInNLEITUNG 


Bei den mechanischen Trennverfabren kann man 
zwei Prinzipien unterscheiden [1] : 


1. das Absetzprinzip 
2. das Siebprinzip. 


Die der Absetzver- 


fahren, bei 


praktische Anwendung 


denen sich feste Teilchen in. einer 


Flissigkeit bewegen, ist vielfach der Erforschung 
Die 
vorliegende Arbeit soll daher zur Klarung der 


der theoretischen Grundlagen vorausgeeilt. 


theoretischen Zusammenhange beitragen. 
Bewegt sich ein Teilchen in einer Flissigkeit, 

Widerstand, d.h. 

Teilchen wird eine Kraft ausgeiibt, die entgegen- 


so erfahrt es einen auf das 
gesetzt der Bewegungsrichtung an ihm angreift. 
Die Bestimmung der Widerstande von Koérpern, 
Flissigkeiten und die 
Ableitung von Formeln fiir die Sinkgeschwindig- 
keiten dieser Kérper in Flissigkeiten beschiftigte 


die Forschung seit langer Zeit. Zunichst wurden 


die sich in bewegen, 


Beziehungen zur Berechnung der Widerstande, 
die bewegte Kugeln in einer Flissigkeit erfahren, 
gefunden. In der Technik hat man es jedoch 


meist mit unregelmissig geformten Kérpern zu 
tun, so dass die abgeleiteten Formeln nur ange- 
naherte Gultigkeit haben. 

Wahrend sich eine frihere Arbeit von Lupwic 
[2] mit der Sinkgeschwindigkeit von Scheiben und 
Zylindern beschaftigte, untersucht die vorhegende 
Arbeit fiir isometrische (regulare) Teilchen, wie 
Wirfel, Oktaeder und Tetraeder, den Zusammen- 
hang zwischen dem Widerstand, den ein solches 
Teilchen in emer Flissigkeit erfahrt, bzw. dessen 
Sinkgeschwindigkeit und der Teilchenform. 


1.1 Mathematische Grundlagen 

Der Widerstand W, den ein in einer Flissigkeit 
bewegtes Teilchen erfahrt, wird allgemein fol- 
gendermassen formuliert [3, 4] : 


: v2 
W = pF pp (1) 


Nach der Abhangigkeit des Widerstandskoeflizien- 
ten & von der Reynoldsschen Zahl Re unter- 
scheidet man drei Bereiche : 

(a) Stokesscher Bereich (0 << Re < 0,2). In 
diesem Bereich gilt ¢ = 24 Re. Die Stokessche 
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Formel fiir den Widerstand von 


(F = dat Re 


durch Einsetzen dieser Werte in Gleichung (1): 


Kugeln 


Ugd v = Ugdpp,y) erhalt man 


W — 8 dv, (2) 


(b) Allenscher Bereich (30 < Re < 300). Hier 
ist wh 10/4 Re. so dass die Allensche Formel fiir 
den Widerstand von Kugeln lautet 

" —— d8 pp va? (3) 
ty Re 


Bereich (Re > 1000), In 


% unabhangig von Re, d.h. 


(c) Newtonscher 
diesem Bereich ist 
wu — const. Damit ergibt sich die New- 


tonsche Formel fiir den Widerstand von Kugeln 


Cy . 


(4) 


Die Gebiete zwischen diesen Bereichen werden 


von keiner Formel exakt erfasst, so dass man auf 


Aus 
Masse 


emes 


graphische Darstellungen angewiesen ist. 
der dynamischen Grundgleichung Kraft 
folet fur Sinken 
Flissigkeit (4 — Auftrieb, 
Widerstand) 


Beschleunigung das 


Teilchens in’ einer 


S Schwerkraft, W 
m 4 S-A-W (5) 
dt 


Hat das Teilchen seine maximale Sinkgeschwin- 
digkeit erreicht, so ist dv ‘dt = 0, d.h. 
We S-A (6) 


Durch Einsetzen der Werte fur eine Kugel folgt 
daraus fiir deren Sinkgeschwindigkeit allgemein 


pet 42 pr — Pry 
x 3 ou Pr 


(7) 


Somit ergibt sich fir die Sinkgeschwindigkeit v, 
einer Kugel 


(a) im Stokesschen Bereich mit ¢ = 24 Re 


d? 


UK (pr - 


187 Pr) fs (8) 


(b) im Allenschen Bereich mit 


*/[ 4g* (pr — pr)? 
ler » ps Os 1 P 
"KN [225 pp ” 


10 \ Re 


(c) im Newtonschen Bereich mit % — const. 


Cy 


(10) 


(11) 


bzw. mit 


1.2 Einfluss des Wandeffektes auf den Widerstand 
und die Sinkgeschwindigkeit 

Die Stokessche Formel (Gleichung 2) fiir den 
Widerstand, Kugel in emer 
Fliissigkeit erfahrt, setzt voraus, dass im Giiltig- 
(0 < Re < 0,2) das 


ausgedehnt ist. 


den eine bewegte 
keitsbereich dieser Formel 
Fliissigkeitsvolumen unendlich 
Bei endlichen Gefissdimensionen muss zur Elimi- 
nierung des Einflusses der Gefasswande auf den 
Stokessche 
erganzt werden. Von den zahlreichen vorgeschla- 
Korrekturfaktoren 


des Wandeffektes wurde hier der von LADENBURG 


Widerstand die Widerstandsforme!l 


yenen zur Bericksichtigung 


[5] verwendet. Danach betragt der Widerstand 
einer Kugel vom Durchmesser d, die langs der 
Achse eines Zylinders mit dem Durchmesser D 
fallt : 
Ww Bandi, (1 2,.4d/D) (13) 
1.3 Definition der Teilchengrésse und Charakteri- 
sierung der Teilchenform 
Die 
Teilchen, die 


Definition der Teilchengrésse ist bei 
Kugelform 


So herrscht 


von der abweichen, 


nicht einheitlich. in der Literatur 
einige Verwirrung beziiglich dieser Definition und 
es besteht eine Vielfalt von Begriffen. 

Arbeit wird zur Festlegung der 
Teilchengrésse nur der von Wape ut [6] vorge- 


schlagene “* wahre Nenndurchmesser 


In dieser 


dy ver- 
wendet. Der “ wahre Nenndurchmesser ” dx ist 
detiniert als der Durchmesser einer Kugel, die bei 
gleicher Dichte dasselbe Volumen hat wie das 
betrachtete, nicht kugelférmige Teilchen. 

Zur Charakterisierung der Teilchenform wird 
hier der ebenfalls von Wave. [7] eingefiihrte 
* Kugelwert ” Ky benutzt; dieser ist definiert als 
das Verhaltnis der Oberfliche Ox einer Kugel zu 
der Oberflache Oy, eines nicht kugelférmigen 
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Teilehens, wobei Kugel und nicht kugelférmiges 
Teilehen dasselbe Volumen, also denselben wahren 
Nenndurchmesser dx haben (Abb, 1a) : 


Ox 


Ky (14) 


OnK 


. oP". 





Wur fe! 


Volumen Ves 


Kugel 
Volumen V, =Vey 
Oberfiladche 0, 
Durchmesser oc =4d, 


a 
Wur fe! 
Volumen Vw 
einbdeschriebene Kuge! 
Durchmesser ¢ 


Kuge! 
Volumen Vi = Vinvw 
Durchmesser c= Cc, 


Ans. 1. Definition des Kugelwertes Ay am 
des Wiirfels 
(a) allgemeine Definition 


Beispiel 


(b) Definition fiir isometrische Teilchen 


Fiir isometrische Teilchen lasst sich noch eine 
andere Formulierung des Kugelwertes Ay an- 
geben. Fir solche Teilchen gilt allgemein [8]: 


(15) 


(d* — Durchmesser der dem isometrischen Teilchen 
einbeschriebenen Kugel). 
Daraus folgt : 
Ox 6 ‘ 
_X — — (fir Kugel), 
J K dx 
Onx 


5 , : , 
; (fur isometrische Teilchen) 
VK 


( 
Qs (16) 
Nach der Definition 


Ve = Vy sein, dh. 


des Kugelwertes muss 


Ox dx = Onk d* 


Aus Gleichung (17) ergibt sich mit Gleichung (14) 
fiir den Kugelwert isometrischer Teilchen 


d* 


Kw dx (18) 


Der Kugelwert isometrischer Teilchen  lisst 
sich also ausdriicken durch das Verhaltnis des 
Durchmessers d* der dem Teilchen ecinbeschrie- 
Kugel 
Teilchen volumengleichen Kugel (Abb. 1b). 


benen zum Durchmesser dg der dem 


2 EXPpeERIMENTELLER TEIL 
2.1 Versuchsanordnung 


Die Versuchsanordnung (Abb. 2) besteht aus 
einem Glaszylinder a mit einem inneren Durch- 
messer D — 76 mm und einer Lange von 1150 mm; 
dieser Glaszylinder dient als Raum fir 


Sinkversuche. Er gegen einen 


die 


kann auch 








App. 2. Schema der Versuchsanordnung. 


Glaszylinder mit einem inneren Durchmesser 
D = 42,5mm und einer Linge von 1150 mm 
ausgewechselt werden. Glaszylinder (a) befindet 
sich auf einer Korkunterlage in einem zweiten 
Glaszylinder (b) mit einem inneren Durchmesser 
D = 100mm und einer Linge von 1150 mm. 
Ein Korkring verschliesst den Zylinder (b) und 
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dient gleichzeitig als Halterung fiir (a). Der Ring- 
raum, der durch die beiden Zylinder gebildet 
wird, wird durchflossen, dessen 
Temperatur durch den Thermostaten (c) auf + 0, 
1°C konstant gehalten wird. Die Sinkstrecke im 
Glaszylinder (a) ist durch zwei Marken festgelegt, 
die (500 + 0,5) mm _ voneinander entfernt sind. 
Die obere Marke befindet sich 450 mm unterhalb 


des 


von Wasser 


untere 
Marke 150 mm vom Boden entfernt ist, um einen 
Einfluss des Bodens auf die Sinkgeschwindigkeit 
auszuschalten. 


Flissigkeitsspiegels, wahrend die 


Bei den Sinkversuchen wurden Tetraeder, 
Wiirfel und Oktaeder aus Kunststoff (p, 1,213 
1.314 ¢ em*) Kugeln Glas 
(pr — 2.930 g em) und Stahl (p7 — 7.760 g cm*) 
in Abmessungen von 0, 95 bis 45 mm verwendet. 


Versuchsteilchen und Versuchsflissigkeiten 


bis sowie aus 


Es wurden die Sinkgeschwindigkeiten einzelner 


Teilchen gemessen. Teilchen, an denen Laft- 








blischen hafteten und solche, die entlang der 


Wandung des Glaszylinders sanken, wurden 
bei der Auswertung nicht beriicksichtigt. Als 
Versuchsfliissigkeiten dienten destilliertes Wasser, 
wassrige Calciumchlorid-, Zucker- und Glyzerin- 
lisung, ferner Mineralél. Die Viskositaéten dieser 
Fliissigkeiten konnten noch durch Anderung der 
Temperatur variiert werden; der Viskositats- 
bereich erstreckte sich so von » 6,4 10-° bis 
7,48 g/cm sec. Die Dichten der Flissigkeiten 


lagen zwischen pp — 0.9087 und 1,314 g ‘cm*. 


3 VERSUCHSERGEBNISSE 

Aus den gemessenen Sinkgeschwindigkeiten 

im Stokesschen Bereich (0 < Re < 0, 2) wur- 
den diese zur Beriicksichtigung des Wandeffektes 
nach der Ladenburgschen Formel korrigiert (siehe 
Abschnitt 1.2) 
durch Einsetzen des wahren Nenndurchmessers 
dx an Stelle von d der Widerstandskoeffizient px 


lasst sich mittels Gleichung (7) 


berechnen, wobei der Index AK ausdriicken soll, 














| 














10 
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4 68yy 2 


a 6 8105 2 


Re, 


Ans. 3. 


Zahl Re, (a Tetraeder, b 


Widerstandskoeflizient ¥, als Funktion der Reynoldsschen 
Wiirfel, c 


Oktaeder, d = Kugel). 


95 





K. ¢ 


dass der Berechnung der wahre Nenndurclimesser 
dx zugrunde gelegt wurde. Auch die entsprechen- 
den Reynoldsschen Zahlen Rex wurden mit Hilfe 
des wahren Nenndurchmessers dx berechnet. 

Fir den Bereich Reynoldsscher Zahlen von 
etwa 0,003 bis 20,000 sind die Widerstandskoefli- 
zienten vd, als Funktion von Reg im doppelt- 
logarithmischen Netz in Abb. 3 dargestellt. Die 
Kurven dx = f(Re,x) fir Tetraeder (a, Kugelwert 
Ky 0,670), Wirfel (b, Kugelwert Ay 0.806), 
Oktaeder (c, Kugelwert Ay 0.846) und Kugeln 
(d, Kugelwert Ay 1) liegen mit abnehmendem 
Kugelwert d.h. je 
Kugelwert ist, desto niedriger liegen die Kurven, 


tibereinander, héher der 
desto niedriger ist auch der Widerstandskoetlizient 
Wy. Lediglich im Bereich von etwa Re, — 20 
bis Rex 140 liegt die Kurve fiir Oktaeder (c) 
iiber der fir Wirfel (b), obwohl Oktaeder einen 
héheren Kugelwert haben. 


3.1. Berechnung der Sinkgeschwindigkeiten nicht 
kugelform ige r Te ilche n Stoke ssche ” Be reich 
(0 < Re < 0, 2) 
Die Stokesschen 
Sinkgeschwindigkeit, Gleichung (8), gilt nur fir 
Sink- 
geschwindigkeit nicht kugelférmiger Teilchen nach 


am 
Formel zur Berechnung der 


kugelférmige Teilchen. Um auch die 


dieser Gleichung berechnen zu kénnen, fuhren 


wir einen Formkorrekturfaktor ein durch die 
Beziehung 

Unq = Avy (19) 
Darin ist vy, die Sinkgeschwindigkeit eines nicht 
kugelférmigen Teilchens, dessen wahrer Nenn- 
durchmesser dy gleich dem Durchmesser d einer 


volumengleichen Kugel ist, die unter denselben 


Tabelle 1. 


des Versuchszylinders, Ky 


Formkorrekturfaktoren K fiir Tetraeder, Wirfel, Oktaeder und Kugeln (D 


R. Cnowpuury und W. Frrrz 


Bedingungen mit der Geschwindigkeit vg sinkt. 

Aus Gleichung (19) erhalt man durch Einsetzen 

von Gleichung (8) und mit d = dx, nach einfacher 

Umformung fiir den Formkorrekturfaktor 
IS » Var 

) UNK (20) 


dx g (pr — pr) 


Mittels Gleichung (20) lassen sich die Formkorrek- 
turfaktoren aus den experimentell ermittelten 
Sinkgeschwindigkeiten Zur Elimi- 
nierung des Wandeffektes, der vom Verhaltnis 
dy D abhangt (siche Abschnitt 1.2), wurden die 
nach Gleichung (20) berechneten Formkorrektur- 


berechnen. 


faktoren jeweils fiir eine bestimmte Teilchenform 
graphisch tiber dy aufgetragen. Die Formkorrek- 
turfaktoren sind linear abhangig, 
dh. K b. Die Konstanten a und b der 


Geradengleichung wurden nach der Methode der 


von dx 


ad x 


Summe der kleinsten Fehlerquadrate von Gauss 
[9] bestimmt. Fir d, — 0 erhalt man den reinen 
Formkorrekturfaktor (— Ordinatenabschnitt 5), 
Die Extrapolation auf den wahren Nenndurch- 
0 gleichbedeutend mit = der 
Stokesschen 


unendlich ausgedehnten Flissigkeitsvolumens. 


messer dx ist 


Voraussetzung des Gesetzes eines 


Zum Vergleich wurden die Sinkversuche sowohl 
in einem Glaszylinder mit 42, 5mm Innendurch- 


messer als auch in einem solchen mit einem 


76,0mm_ durehgefihrt. 
Messergebnissen beider Versuchs- 


Innendurchmesser von 
Die 


reihen ermittelten Formkorrekturfaktoren zeigten 


aus den 


gute Ubereinstimmung (Tabelle 1). 
Zwischen dem Kugelwert Ay und dem Korrek- 
turfaktor A besteht eine einfache mathematische 


Beziehung. Wird der Logarithmus des Kugel- 


wertes als Funktion des Formkorrekturfaktors 
Tnnendurchmesser 


Kugelwert) 





Tetraeder 


Ky 


V ersuchsrethe 
o.670 


1(D 42.5 mm) O85 


2(D 76.0 mm) O85 


Mittelwerte O85 


Wirfel 
Ky 


Oktaeder 
AK u O.846 


Kugel 
O.806 K Ww 1 


O91 O04 


O05 O04 


O.v2 O04 








Sinkversuche mit isometrischen Teilchen in Fliissigkeiten 


d.h. 
Konstanten 
Tabelle 1: 


eine Gerade, 


Die 
der 


sich 
b (Abb. 4). 
den Daten 


aufgetragen, so ergibt 
log K " ak 


errechnen sich aus 


Mw O7 
08 
A 


0,6 





05 1 
06 07 


69 140 


Ans. 4. Zusammenhang zwischen Kugelwert 


Formkorrekturfaktor A. 


Ky und 


Damit 
durch 


1.159, bd 
1159 K 
folgt : 


1,159. 


woraus 


wird log Ky 
1,159, Umformung 


K 1 0.862 « log Ky (21) 


Nach dieser Gleichung, die fiir Tetraeder, Wiirfel. 
Oktaeder und Kugeln mit einer Genauigkeit von 

1% erfillt ist, lassen sich aus den Kugelwerten 
die Formkorrekturfaktoren und damit die Sink- 
geschwindigkeiten vy, dieser Teilchen im Stokes- 
Bierech nach der Gleichung 


schen berechnen 


(22 


die sich aus Gleichung (19) durch Einsetzen von 


Gleichung (8) mit d d,. ergibt. 


3.2 
kugelformiger Teilchen im Newtonschen 
(Re > 1000) 


Die Konstante Ay der Newtonschen Forme! 
fir die Sinkgeschwindigkeit, Gleichung (12), kann 


Berechnung der Sinkgeschwindigkeiten nicht 


Bereich 


unmittelbar aus den experimentell ermittelten 
Sinkgeschwindigkeiten berechnet werden, wenn 
in Gleichung (12) fir d =—d, und fiir vy die 


Tabelle 2. 


Sinkgeschwindigkeit vy, der nicht kugelfOrmigen 
Es zeigte sich, dass fiir 
Ky 


vom wahren Nenndurchmesser dx, ist. 


Teilchen gesetzt wird. 


eine bestimmte Teilchenform unabhangig 

Die aus den Sinkgeschwindigkeiten berechneten 
Werte von Ky sind in Tabelle 2 aufgefiihrt und in 
Abb. 5 als Funktion des Kugelwertes Ay, darge- 


stellt. 








06 i 
0,6 0,7 0,8 0,9 


Kw 





1,0 i1 


Anns. 5. Konstante Ky der Gleichung (12) als Funktion 


des Kuygelwertes Ky- 


Die Kurve der Abb. 5 hat in erster Naherung 
Verlauf ; kann Parabel 
dritten Grades angesehen werden mit der allge- 
bx + ca* + d2*. Die 
Berechnung der Konstanten dieser Gleichung aus 
den Werten der Tabelle 2 ergibt fir Ky: 


Ky 


parabolischen sie als 


meinen Gleichung y = a 


oo 


2.22 —8,04Ky —l11Ky+3,71Ky (28) 


Mit Hilfe dieser Beziehung, die fiir Tetraeder, 
Wirfel, Oktaeder und Kugeln mit einer Genauig- 
keit von 1% erfillt ist, 
Kugelwerten Ay die Konstanten Ay fir diese 


lassen sich aus den 


Teilehen berechnen. Deren Sinkgeschwindigkeit 


ergibt sich aus folgender Gleichung 


—_ Ky Je PT ~ Pr dx 


Konstante Ky der Gleichung (12) fiir Tetraeder, Wiirfel, Oktaeder und Kugel 





Teilchenform Tetraeder 


Kugelwert Ay 0.670 


Ky 0.79 OO 0.08 


Wiirfel 


0.806 


Oktaeder Kugel 


0,846 1 


+ OO1 1,09 + O01 1,79 + O01 
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Ky 
Ky 


die aus Gleichung (12) durch Setzen von vx 

und d—d, folgt. Im Bereich zwischen 
Stokesschen und Bereich 
lassen sich keine Beziehungen ber den Zusammen- 


UNK 
dem 
dem Newtonschen 
hang von Sinkgeschwindigkeit und Teilchenform, 


die durch den Kugelwert charakterisiert ist, 


aufstellen. 


Erklarung der Formelzeichen 


(Die Formelzeichen bedeuten Groéssen. Die bevor- 


zugten Einheiten sind in Klammern angegeben) 


Cy KRonstante der Newtonschen Widerstandformel 


D Durchmesser des Messrohres om 
d Durchmesser ciner Kugel em 
d* Durchmesser der einem tsometrischen 

Teilchen einbeschriebenen Kugel 
aK wahrer Nenndurchmesser Durchmesser 
emer Kugel, die bei gleicher Dichte das- 
selbe Volumen hat wie em nicht khugelfor 
thiges Teilchen 
ingest romte Querschnittflache emes 
Pealchens en? 
birdbeschleunwuny em sec? 


Formkorrekturfaktor 
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Vapour-liquid equilibria at reduced pressures : n-butanol—ethyl benzene ; 
aniline=p-cymene 


S. R. M. Exviis* and M. Razavieour 
( Received 19 February 1959) 


Abstract The increasing importance of vacuum distillation in the petroleum, chemical and allied 
industries has drawn attention to the need for accurate vapour liquid equilibria results at 
reduced pressure. 

Vapour liquid equilibria results are presented for the azeotropic systems n-butanol ethyl 
benzene, and p-cymene aniline at pressures from 50 to 760 mm Hg. \ test is made for the 
thermodynamic consistency of the data, and the limitations of the correlating equations for 


isobaric conditions are discussed. 


Résumé—L ‘importance croissante de la distillation sous vide dans lindustrie du pétrole, Vindustrie 
chimique et autres industries connexés a attiré lattention sur le besoin de résultats précis des 
équilibres liquide- vapeur a pression réduite. 


les résultats des équilibres liquide-vapeur poru les systémes 


Les auteurs présentent 
azéotropiques n-butanol — éthylbenzéne et p-isopropyloluéne 
50 a 760mm Hy. Hs veérifient la consistance thermodynamique des données et discutent des 


limites des équations correspondantes pour des conditions isobares. 


aniline pour des pressions de 


Zusammenfassung_ -Die steigende Bedeutung der Vakuumdestillation in den Erdol-, chemischen 
und verwandten Industrien erfordert genaue Kenntnis der Dampflliissiggleichgewichte bei kleinen 
Drucken. 

Dampfiliissiggleichgewichte werden fiir die azeotropischen Systeme n-Butanol — Athylbenzol 
und p-Cymolanilin bei Driicken von 50 bis 760 Torr mitgeteilt. Die thermodynamische Konsistenz 
der Daten wurde nachgepriift und die Begrenzung der Beziehungsgleichungen fiir isobare 


Bedingungen diskutiert. 


APPARATUS 
Tue still used was the recirculating type, pre- 
viously described [1], and the pressures ranged 
from 100 to 760 mm Hg. 

When the still was operated at 50 mm it was 
necessary to modify the heater unit for even 
boiling. The clearance between the inner and 
outer walls of the bayonet heater unit was re- 
duced from 4mm to Imm. The tube joining the 
top of the boiler to the Cottrel coil was reduced 
from 8mm to 6mm. These modifications had 
two taain advantages ; firstly, the velocity of the 
liquid over the heating surface was increased, thus 
improving the rate of heat transfer; secondly, 
when any “ bumping” occurred, the volume of 
liquid forced from the boiler into the still was 
smaller, and had a comparatively small effect on 
the smooth running of the apparatus. 


PressurE CONTROL 


When the still was operated under vacuum the 


pressure was controlled by the control system 


shown on Fic. 1. The actual automatic controller 
uses a sintered disk glass and was similar to that 
suggested by M. L. Freperick [2]. The chambers 
D and B were two cylindrical glasses 3}in. high 
and lin. diameter, connected by a U-shape tube, 
H, }in. diameter. The chamber D could be closed 
off from the pump by the stopcock E while the 
stopcock C disconnected the still from the pump. 
The required pressure was obtained in chamber D 
by leaving open the stopcocks C and E until the 
manometer, M, indicated the required pressure 
and then stopeocks C and E were closed. The 
pressure in chamber D was the operating pressure 
of the still and an increase in pressure in the still 
resulted in a movement of mercury in the manostat 
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TO STILL 


MANOME TER 
TO ATMOSPHERE o<f MANUALLY 
VALVE 

T VACUUM PUMP 


- 
Ls 3) 


Fie. 1 


tube so as to open the sealed disk, Q, between B 
and F. When the pressure in the chamber D was 
equal to that in the still, the mercury rose thus 
sealing the disk again. 

\t the end of the run, cocks C and FE were open 
that 
An auxiliary manometer, 


ed in order to bleed air into the still. so 
samples could be taken. 
S, which indicated the pressure of chamber B was 
installed to prevent disturbing the apparatus. 
Before opening the cocks, air was let through by 
opening the valve Z until the pressure in B became 
nearly equal to that in D, the cocks C and E wer 
then both opened together. A further advantage 
in using the manometer S, was that the pressure 
in chamber D could be altered during the runs if 
required, 

Marvrertan Usen Mi 


AND rHoDS Oo} 


ANALYSIS 
Analar grade ethylbenzene and n-butanol were 


purified in a 30 plate Oldershaw column at 


Pressure 


PRESSURE REGULATOR 


WERATED 


control system. 


atmospheric pressure while aniline and p-cymene 
purified the 
apparatus at The physical 


were by fractionation in same 


reduced pressure. 
properties of the purified materials are tabulated 
in Table 1. 

The 


p-cymene 


vapour pressures of ethylbenzene and 
calculated using the 
equation reported by the National Bureau of 
Standards [6]. Data by Sru.s. [7] were used for 
n-butanol; the log of the 


were Antoine 


Vapour pressure Was 
plotted against the log of the vapour pressure of 
pure water at the same temperature. Equations 
were then derived from the resulting straight line 
plots. Vapour pressure results for aniline were 
taken from TimMerMaANs [3]. 

For the analysis of the binary mixtures the 
refractive index-composition plots at 20°C were 
ethyl 

mixtures, 


benzene n-butanol and 


The 
An Abbe refractometer was 


determined for 


aniline. p-cymene results are 
presented in Table 2. 


used for refractive index measurements. 


Table 1. 





Refractive indes 
uf 20° 
Materials . 


Purified Literature 
1-4958 
11-3993 
1-5862 
14903 


Ethylbenzene $058 [3] 


n-Butanol 3993 
Aniline 5862 


p-cymene 4908 


Purified 


O-8671 
O-8097 
10219 
O-8571 


Boiling t 
Density 20.4 0 = poin 


Literature Purified 


136-1 
117-7 


184-1 


Literature 


O-8671 [3] 
0-8097 [3] 
1-0218 [3] 
O-8573 [6] 177-7 


136-1 [3] 
117-7 [3} 
184-35 [3] [4] 
177-2 [6] 








Vapour liquid equilibria at reduced pressures : n-butanol_ethyl benzene ; aniline p-cymene 


Table 2 butanol on Fic. 2 shows that a minimum boiling 





azeotrope exists at all pressures investigated, with 
p-cymene aniline Ethyl benzene p-cymene the azeotrope composition becoming less rich in 
n-butanol as the pressure is reduced. It will also 


Mol. ° + Mol. ° be noted that the concentration of butanol in 


o 


= RA. at 24 Rl. at 2 
preymene n-butanol the vapour decreases as the pressure is reduced. 


Fig. 3 shows that there is a minimum boiling 
Oa ost. 58 


SIS = Wade bate 
18-88 558: 1876 concentration of p-cymene in the vapour increases 


azeotrope in the system aniline—p-cymene. The 


28-10 ! 1800 as the pressure is increased. 

Ste Ot) 53: ; t 4751 

48°45 ° ‘ wis 

5850 ey We f $528 

6s 0 Pr jot $455 
Tuam) ies. tS 35 bs 

SANG 50 ut 2233 

1a bem sf bliss 
loo BSUS) 





EXPERIMENTAL RESULTS 


Vapour liquid equilibrium results for the system 
ethyl benzene-v-butanol for 760, 500, 300, 100 
and 50mm Hg are presented graphically on 
Fic. 2 while results for the system aniline—p- 


cymene are shown on Fic. 3. Typical experimental 
results are shown in Tables 3 — 6. 


The equilibrium results for ethyl benzene—n- MOL YMENE IN LIQUID « 


Vapour liquid compositions for aniline p- 


cyvinene, 


AAPOR 


CORRELATION OF EQUILIBRIUM 
ReEst LTs 


N 


Activity coeflicients were calculated by the 


BUTANCK 


equation 


° 
lef 


7Yy 


Ii» 
Po 


moe 


where 7 = total pressure 


P, = vapour pressure of pure component at 
boiling point. 
icitennaieescilatiiiantDatiaDemdieas mole fraction of component in the 
0 46 50 60 © 8 GO Wo I 
etait Oc tiie en cite vapour phase. 


Vapour-liquid compositions for n-butanol mole fraction of component in the 


ethylbenzene. liquid phase. 
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Table 3. Ethylbenzene (e)-n-butanol (b) 





a Composition Activity coefficients 
Total 


pre ssure 


Boiling 
point Liquid Vapour 


(mm) eC) % log Y> log Ye 


117 WS 05-55 O-0003 O-3973 
116-7: we 0-00 O-OOLO O-BTST 
116-2 ' S405 O-O071 O-3482 
115-8: . 80-55 OO1s4 Osli2 
115-85 Sot 7600 O-0207 0-2783 
113-85 2s! T2585 O-0r289 2502 
158 ities TOS O-0506 0-2068 
1160 itt 66-80 O-0580 1882 
116-3. 57-3: O70 OOTz O-1619 
116 $4 62-10 O-OSTO 00-1437 
116-0” i Goon o1tie lis2 
117-50 2-4 57-05 O- 1326 O-0966 
118-15 : 5500) 00-1630 O-OT3T 
119-05 B1- 51-00 O-1918 O-0577 
120-30 2: 8-05 0-2308 0-038 1 
121-80 0-7. SS O-2sSel O-O232 
125-00 . 38-30 O-3147 0-0086 
126-80 20-05 Oss20 OWIT2 
i20-75 56 22s O3BTST O-O0r28 
131-45 f 15-80 O-3736 O-0003 





Table 4. Ethyl benzene (e)-n-butanol (b) 





Total Composition Activity coefficients 
ola 
Run Boiling 
pressure 
no. point Liquid apeur 
(mm) log yy, log y, 
(¢) a» Mp 


100 67-95 96-00 W005 O-O188 0-5042 
100 67-35 93-25 84-40 WOLS2 O-4825 
100 66-35 OSS TO-00 0-282 O-4527 
100 65-85 5M) 73-00 00-0309 0-3946 
100 65-00 “) 68-00) 0-0400 O-3736 
100 64-55 5-30 GS-85 O-O515 O-3B351 
100 64-20 Ww O35 00-0658 0-2989 
100 64-00 65 56-50 O-0837 O-2581 
100 63-65 10 54-50 O10 2324 
TOO 63-80 10 52-15 O-1234 O-2012 
00-7 O3-60 Mw 40-05 | 00-1553 00-1730 
100-3 63-80 270 47-75 0-1880 O- 1448 
100 Os-70 BO. 45-05 O-2380 00-1166 
loo O4-35 B2-6 1-00 2541 0-054 
Oo O445 ° 41-00 2097 00-0803 
100 64-25 . 39-25 0-3958 OOT18 
100 65-40 . 30-35 04171 0-429 
VO-S5 65-85 2: 82-15 00-4899 0-0332 
TOO 67-80 . 25-60 0-6166 00-0245 
100-6 69-90 . 16-65 0-6265 O-O147 
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Table 5. 


p-cymene (c)— Aniline (a) 





Total 


| pressure 
i 
(mm) 


TOO 
760 
760 
760 
760 
760 
760 
760 
TOO 
Ti 
7h 
THO 
TH 
760 
TOO 
THO 
760 
7TH 
760 
TOO 
THO 


Boiling 
point 
(C) 


174-80 
174-05 
173-60 
173-20 

73-00 
172-90 
172-85 
172-80 
172-85 
172-85 
1730 
173-10 
173-75 
174-30 
174-90 
175-55 
176-10 
177-30 
178-75 
179-20 
180-60 


Composition 


Liquid 
He Y, 


so-00 
83-00 


93-00 
87-50 
83-00 79-25 
79-50 76-00 
73-25 70-80 
Os-O0 66-80 
64-00 63-00 
0-50 

50-50 

58°25 

5240) 

W530 

450 

OO 

36-50 

32-50 

28:75 

22-40 

16-40 

14-40 

SSO) 


V vapour 


Activity coefficients 


log y, log Ye 


O-O1l2 
00-0163 
00-0193 
00-0242 
0-03 16 
0-040) 
O-O476 
O-O4AS6 
O-0538 
O-0554 
O-OS44 
O-0945 
1250 
1465 
1774 
O-1919 
2217 
0-2604 
O-P297 
0-3300 
OSS 


00-2691 

0-2288 
00-2109 
O-1972 
O-1702 
00-1486 
O-1342 
0-1300 
00-1250 
(12382 
01-0930 
OS19 
00-0641 

O-0552 
O-O406 
O-O3T0 
00-0320 
O-O270 
0-0006 
OO“ 
OLS! 





Table 6. 


p-cymene (c) — Aniline 


(a) 





Run 
no. 


Total 
pressure 
(mm) 


99-5 
100-0 
100-0 
100-4 
100-4 

90-8 

99-5 

99-25 

99-6 

99-9 
100-3 
100 
100 
100-3 
100-1 

99-7 
100 

99-75 

99-5 
100-3 





Boiling 
point 
(©) 


107-30 
107-00 
106-80 
106-70 
106-80 
106-40 
106-30 


109-00 
109-15 
109-75 
110-70 
111-40 
112-85 
114-35 


Compositions 
Liquid Vapour 
ae Ye 
93-00 
89-50 
85-40 
80-50 
76-00 
69-10 
65-30 
60-00 
55-75 
50-70 
46-80 
42-00 
38-50 
29-20 
24-90 
20-80 
16-75 
13-40 
10-40 

6-60 


SO-O0 
s6-10 


20-00 


Ietivity coefficients 


log Ye log Ya 


O-O135 
OOoLT2 
O-O198 
O-256 
0-0306 
O-O452 
00-0566 
0-0668 
0-O847 
00-1009 
01143 
0-1234 
O-1519 
O-2212 
O-2516 
(2732 
00-3167 
O-3507 
0-3676 
O-MOAT 


o-s602 
O-3444 
O-3147 
0-287 1 
O-2525 
02248 
02054 
01842 
0-1639 
0O-1401 
O-1190 
0-1025 
O-O875 
00-0529 
O-O62Z0 
0-0528 
0-0384 
0-0354 
0-0242 
0-282 
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== 


Fic. 4. Activity coefficient curves for n-butanol-ethyl- 


benzene, 


Activity coefficient curves for n-butanol ethyl 
benzene at 760 mm Hg. 


~~ 
a 


a ae ee ee 
MOLE 9, 9- BUTANOL IN LIQUID 
6. Activity coefficient curves for n-butanol ethyl- 
benzene at 100 mm Hg. 


MOLE % 2-BUTANOL 


Repuice and Kisrer plots for n-butanol ethyt- 
benzene, 


Ethylbenze ne n-butanol 

Fig. 4 shows plots of log activity coeflicient 
against concentration for different pressures, It 
can be seen that for a given concentration, the log 
of the activity coefficient is increased as the 
pressure is reduced. Figs. 5 and 6 show a com- 
parison between experimental and van Laar plots 
for 760 and 100 mm pressure. 

On Fig. 7 of the log of the ratio of the activity 
coefficients is plotted against composition for 
760 mm and 100 mm pressure. Similar plots were 
made for other pressures. In Table 7 a compar- 
ison is made of the areas above and below the zero 


line. 


Table 7. 





trea above the | Area below the . 
Pressures , ©. Error 
zero line zero line 


10-05 12-82 
11-35 125 
12-85 13-1 
14:55 14-14 
1540 14-4 








Vapour- liquid equilibria at reduced pressures : 


The percentage error in Table 7 is expressed as the 
decrease or increase in area below the zero line 


compared with that above the zero line. 


p-cymene—analine 

Log activity coefficient plots are shown on Fig. 
8, and again it can be seen that the log of the 
activity coefficient increases as the pressure is 
reduced for a given composition. When a com- 
parison is made between experimental and Van 
Laar plots the agreement is good. 


Table 8. 





Pressure trea above the Area below the 


‘> erro 


(mm Hg.) zero line wero line 
o4s 
10-16 
300 11-22 
100 sl 
Mw 115 
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700 
Sow 


O4 

o-oo 
O79 
OSS 


~- 
175 


10-09 
11-08 
1-00 
1b 





Table the 
above and below the zero line when for varying 


S summarizes area measurements 


pressure the log ratio of the activity coeflicients is 


n-butanol ethyl benzene ; aniline p-cymene 


plotted against composition. The deviations from 


log "1 dx, 0 
Y 


are much less for p-cymene—aniline than for 
n-butanol-ethylbenzene. 
Discussion oF RESULTS 


Comparing Tables 7 and 8 it can be seen that 
there is a greater deviation between the experi- 
mental and correlated Repiicn and Kister [8] 
than for 


due to a 


results for ethylbenzene—n-butanol 


aniline-p-cymene. This is mainly 
significant difference in the boiling ranges of the 
two systems, although at 760 mm Hg they have 
the same order of non-ideality. At this pressure, 
the boiling-point difference for ethylbenzene N- 
butanol is 15°65°C compared with 7-8°C for 
aniline-p-cymene. 

It would therefore appear that as the derivation 
of the Redlich and Kister correlation is based on 
the isothermal-isobaric Gibbs-Duhem equation, 
the application of the Redlich and Kister method 
to isobaric systems is limited to low boiling 
ranges, 


Herincron [11] and Tuwssen [12] have related 





MOLE kp -CYMENE IN LIOUID x, 


Fie. 8. 


Activity coeflicient curves for aniline-p-cymene. 
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the variation of activity coeflicient with tempera- 
ture (at constant composition and pressure) to the 
heat of mixing, and introduced a correction factor 
into the Redlich and Kister equation for isobaric 
conditions. 


integral heat of mixing for a mole of 
the mixture 


slope of the boiling-point curve 
boiling-point of the mixture in °K, 
H, dT 


RT? dx 
Dopce [9] as B, “the temperature correction 


is referred to by Ist and 


The term | 


factor for isobaric experimental data.” 
In order to determine B values, heats of mixing 
were calculated by the method of Orumer [10], 
log *” 
Ye 
ethylbenzene-n-butanol system and it was found 


B| vs. % Was then plotted for the 


that at 760 mm Hg the percentage error was re- 
duced from 18-4 per cent to — 6-4 per cent, while 
at 100mm Hg it was slightly increased from 

3-0 per cent to 41 per cent. 

Since the calculated heats of mixing are not very 
accurate, as no information was available for the 
variation of the latent heat of vaporization with 
temperature, an attempt was made to use calori- 
metrically determined heats of mixing by Brown 
Fock [13] for the 
butanol at 45-00°C. The assumption was made 


and system benzene n- 

that the heats of mixing for ethyl benzene—n- 

butanol are the same as the experimental values 

for benzene - n Fig. 9 gives a plot 

of oe Ye a| vs. a 
Ye 

n-butanol with a percentage error of only — 0-9 per 


cent at 760mm Hg. and — 50 per cent at 
100 mm Hg. As the boiling range is only 6-3°C at 
100 mm Hg the uncorrected Redlich and Kister 
equation should be applicable. 


butanol. 


for ethylbenzene 


O8, 


0 3 3B 46 5 6 DH BH BD Ho 
MOLE % »-BUTANOL 
Plot of equation (1) using experimental heats of 
mixing. 


It would therefore appear that the limitation 
of low boiling ranges imposed when applying the 
Redlich and Kister equation to isobaric data, 
may be overcome to a large extent by the use of 
equation (1), provided that accurate heat of 
mixing data are available. 


NOTATION 
« = composition of the lighter component in liquid 
y — composition of the lighter component in vapour 
T — boiling point on “K. 
Po, Po, 
Hyp 
R -~ gas constant 


vapour pressure of pure components 


integral heat of mixing (cal mole) 


y = activity coeflicient 


- total pressure 


Subscripts 
1 = lighter component 
2 = heavy component 
a = aniline 
¢ = p-cymene 
b = n-butanol 


¢ = ethyl benzene 





[10) 
(11) 
[12] 
(13) 


Vapour-liquid equilibria at reduced pressures : n-butanol-ethyl benzene ; aniline—p-cymene 
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Abstract 
tangential flow of an isothermal, incompressible viscous fluid in the annular space between two 


An exact solution of the equations of change is given for the unsteady state, laminar, 


evlinders one or both of which may be rotating. The transient behaviour of the velocity profiles 


is set forth graphically. An approximate solution for short times is also given, These results may 


be of interest in connexion with problems in lubrication and viscometry 


Résumé 


tangentiel, laminaire non stationnaire et isotherme d'un fluide incompressible et visqueux dans 


Les auteurs donnent une solution exacte des equations déchangé pour un écoulement 


espace annulaire compris entre deux cylindres dont un seul ou deux cylindres peuvent étre 


rotatifs. Le comportement transitoire des profils de vitesse est présenté graphiquement. Une 


solution approximative est donnée pour des temps courts. Ces résultats peuvent étre intéressants 
par leur rapport avec les problémes de lubrification et de viscosimetrie 

Zusammenfassung—Fiir die Anderungsyleichung der nichtstationdiren laminaren ‘Tangen- 
tialstrOmung einer isothermen inkompressiblen zihen Flissigkeit in einem Ringspalt zwischen 


zwei Zylindern wird eine exakte Losung angegeben und zwar unter der Annahme, dass der eine 


oder beide Zylinder rotieren. 
graphisch fortgesetzt. 


Kine Naiherunygzlésung fiir kurze Zeiten wird ebenfalls mitgeteilt 


Das nichtstationire Verhalten der Geschwindigkeitsprotile wird 


Die 


Ergebnisse sind fiir Probleme der Schmierung und Zahigkeitsmessung von Interesse. 


l Tue Frow Eeuvartions 


THe equations of continuity and motion for the 
flow of a compressible, isothermal fluid are [1, 3] : 


dp ot 


(v- V)v] 


(V-pv)=0 
Vp —(V.-r) eg 


in which p is the fluid density, v is the local 
velocity, p is the static pressure, g is the external 
body force per unit mass and ¢ is the time. The 
stress tensor rt is related to the velocity gradients 
for Newtonian flow thus [1, 3] : 


T we [Vv + (¥v)t] + en (V- vs (3) 
in which »« is the coeflicient of shear viscosity, 


§ is the unit tensor, and (Vv)f is the transpose of 
(Vv). Special cases of the above equations in 


cylindrical co-ordinates are used as the starting 
point in the following developments. 
The 


implies not only that there is no impressed 


assumption above of isothermal flow 
temperature field, but that in addition the viscous 
dissipation term (+: Vv) in the energy balance 


equation is negligible. 


EXPRESSION FOR 
DISTRIBUTION 


EXact 
VELOCITY 


rHE 


A fluid is contained in the annular region 
between a pair of infinitely long co-axial cylinders 
with radii « R and R (see Fig. 1). At time t<0 
the fluid within the annulus is at rest. At ¢>0 
the inner and outer cylinders are rotated with 
constant angular velocities 2; and 2, respectively. 
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For an incompressible fluid in this situation d/19) (4) 
v v, = 0 and v, is a function of r alone so that dé | E dE 
equations (2) and (3) become [1] : 


r 


with boundary and initial conditions: at 
[22 (rey (4) db xa; at &é=1,¢=—1—a; and at +r 
ro d@ = 0. 
This is to be solved with the boundary and 
initial conditions: at r x R, v, «kK RQ,; at 
r= R,vg = RQ,; att < 0,0, = 0. This problem 
is discussed by Kozeny [5] but his solution is 
incorrect in that it does not satisfy the initial 
and boundary conditions, The analogous problem 
for axial flow in an annulus has been solved by 
MULLER [6]. 
It is convenient to introduce the dimensionless 
quantities : 


radial co-ordinate 
pR*® time 
RA Q — tangential velocity 
as 2/A2 angular velocity (8) L ‘ ; 


ore ) () The » ’ > . 
where A & =<6 » The proble m may then Fig. 1. Annular region in which tangential laminar flow 


be restated as: is occurring. 


The quantities (1 x) b,/n 
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The steady state solution at + = & is well- 


known : 


-~a(l «?) 
2 


¢..(€) ' 3 


l1—« 


LS 
1 
Hence the unsteady state solution has the form : 


Hf. rT) = $.(€) — Hf €, 1) 


the function 6/{£, +) being the transient contribu- 
tion, which must satisfy equation (9) but with the 
at & a, &, 


| €1 = At — BE" (10) 
(11) 


boundary and initial conditions : 
d, = O; at r= 0, 4, = $,,( 6). 

A transient function of the form 4/(€, 7) 
f(€) gr) is assumed. Application of the method 
of separation of variables gives : 


1 dg . ld (ef - 
gdr ede | I) (12) 


hoi BE) Z,(by£) € dE 


1 
| [2,(b,€)}* € dé 
- Ab, [Zo(b,) - 
4 [Z,%(b,) 


The integrals in equation (16) are evaluated by 
use of the formulae in Jaunke and Empe [4]. 
The second expression may be simplified by (i) 
inserting the expressions for A and B from 
equation (10), (ii) using the relation 


Jz) Yj) — Ja) ¥,(z) = — 2 ae, 


and (iii) using the defining equation for the ),. 
Thereby it may be shown that 


Z(b,) = — (2/mb,) [J(b,~) J ,(b,) | 
Z(b,,«) = — (2/mb,) (1/ x) 


(17) 
(18) 


x? Z,(b,, x)] > 
x? Z4b,x)] 


F. Curtiss 


where —b*® is the separation constant. The 
solution of the two ordinary differential equations 
thus generated is straightforward, and one 
obtains : 


of €, r) 
wherein 
Z,(b,,€) = J(b,€) V(b, «) — Jb, «) ¥4(b,€) (14) 


The functions Z,(b,£) automatically satisfy the 
boundary conditions that Z,(b,€) = 0 at & = x; 
the requirement that Z,(b,é)=0 at €=1 
determines the eigenvalues b, (some values of 
b, are given in Table 1). The C,, are found from 
the initial condition : 


= C,, (exp - 


by? 7) Zy(byf) (18) 


2 C,, Z,(6,€) 


¢.,(€) 


Multiplication of equation (15) by Z,(6,,€) & dé 
and integration from « to 1 gives for C,, 


Bb,,-* [Zo(b,,) — Zo(b, *)] 


and 


wd (b,, )[a- — a) J,(b,«) + xed (6, )] (19) 


C(x, 
al Xa) [J,%b,«) — J,%(b,)] 


Because of the linear dependence of C,, 
may write : 


on a, we 


Cn, a) = (1 — a) CL, 0) + «Cw, 1) (20) 


Hence the values of C,, for any « may be obtained 
from the tables for « = 0 (inner cylinder station- 
ary) and « = 1 (outer cylinder stationary) given 
in Table 2. 


The complete expression for the unsteady-state velocity profiles is then : 


HE. =|! 


— a(1 a x?) 
54) e [2 
7 SY Aen) (= a) J,(b,«) + 


Je 


ax JS s(by)) [Jilbn€) ¥ilbax) — Ji(bn) Vilbné)) 
[Ji%(bnx) — J1%(,)] 


b,? 7r (21) 
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——— to establish the steady-state distribution may be 
Profiles for a=O Profiles for a=! illustrated by considering four typical examples 
with R = 10cm: 





(i) for water (u/p = 10~-* m*sec!), ¢ = 10*r (sec) ; 
(ii) for air (u/p = 1-3 « 10-5 m*sec~!), t = 7007 
(sec) ; 
(iii) for a typical lubricating oil 
(u/p 10-4 m*sec™!), t = 100+ (sec) ; 


(iv) for a concentrated aqueous solution of 
glycerine (u/p = 5 « 10-4 m*sec™!), t = 20r. 





When += 1 steady-state has virtually been 
established. 


8. APPROXIMATE EXPRESSION 
FoR Vetociry DIstTRIBUTION 
FoR Suort Times 
Equation (21) for the unsteady-state velocity 
profiles is particularly appropriate for large values 
Fie. 2. An illustration of the behaviour of the unsteady- of ve For short times an alternate solution can 
state velocity profiles for « on be derived. 
The p-multiplied Laplace transform of equation 


These dimensionless velocity profiles are plotted (9) ts: 


in Fig. 2 for « = 0 and « = 1. Velocity profiles dfid 
for other values of « may be obtained by adding po dé é dé 
together the velocity profiles for x — 0, 1 according 
to the same rule obeyed by the C,, in equation (20). 

The order of magnitude of the times required 


(68)| (22) 


which is to be solved for the boundary conditions : 
at € Ko ax p;até 1¢ (1 — a«)/p. 


The solution is: 


{(xx/p) Ky(p**) +((1 — 2)/p) Ky(p! 2p} Dip' 2a) +f — 2) pl Ly pe) + (ex p) Lp *)} Kyp'2e) 


, > ’ Ss + (23 
1p) Fa) Kp *) — 1p *) Ky(p' *«) , 


For large p the hyperbolic Bessel functions may be replaced by the asymptotic expressions 
I(r) = (1 2er)'* exp rand K,(r) = (+ 2r)'* exp — r. The resulting expression can then be expanded 
in a manner analogous to that described by Carstaw and Jarcer [2]. An inverse Laplace 
transformation then leads to : 


<.. om erfe 


(¢)! 24( £, r) _ 2n(1 —«)+(2 -£ -) 4 


2y/r 


. 2n(1 — x) + (1 — €) 


2yr 


2n(1 - x) + (1 + € — 2x) 


2yr 


— erfe 


(24) 


which is valid in the limit as + goes to zero. Equation (24) may be further simplified by 
replacing erfe x by a‘ *r' exp —r*(valid for large x) and using only the first term in the 


summations. 
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NOTATION 
Scalars italic 
bold-face roman 
bold-face Greek 


Vectors 
Tensors 
functions defined in equation (10) 
separation constant 
expansion coefficients 
gravitational acceleration 
Bessel function of the first kind (m' order) 
static pressure in Section 1; variable intro- 
duced in Laplace transform in Section 3. 
radial co-ordinate 


y 
ZnO 
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R 
t 
a 
wv) 
a 
£) 
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radius of outer cylinder 


= time 


local velocity vector 

tangential velocity 

Bessel function of the second kind (m*® order) 
Tl on €) Y¥y(b, «) — Fy(d, «) Yyg(b, €) 


- dimensionless angular velocity (equation 8) 


unit tensor 

ratio of radius of inner cylinder to that of 
outer cylinder 

coefficient of shear viscosity 
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Abstract—It is shown how the pressure distribution for isothermal, tangential flow of an 
incompressible fluid in an annulus may be calculated by first solving the problem for a compressible 
fluid and then taking the limit as the deviation from incompressibility vanishes. An analytical 
expression for pressure vs, the radial co-ordinate is obtained for either or both of the cylinders in 
motion 
Résumé Les auteurs montrent comment calculer la distribution de pression pour lécoulement 
isotherme tangentie! d'un fluide incompressible dans une partie annulaire, en résolvant d’abord 
le probléme pour un fluide compressible et en prenant ensuite la limite quand la déviation due a 
lincompressibilite disparait. Une expression analytique de la pression par rapport a la coordonnée 
radiale est obtenue pour Pun ou pour les deux evlindres en mouvement. 


Zusammenfassung —Es wird gezeigt, wie die Druckverteilung fiir cine isotherme Tangential- 
stromung einer inkompressiblen Flissigkeit in einem Ringspalt berechnet werden kann, indem 
man zundichst das Problem fiir cine kompressible Flissigkeit gelést hat und dann als Grenziiber- 
gang die Abweichungen von der Inkompressibilitit zum Verschwinden bringt. Ein analytischer 
Ausdruck fiir den Druck als Funktion der radialen Koordinate wird sowohl fiir den Fall, dass der 
eine oder dass beide Zylinder in Bewegung sind, erhalten. 


It is desired to compute the pressure distribution 
which prevails in the system studied in Part I 
when the velocity distribution is given by ¢_ (£) 
in equation (10) for an incompressible fluid. The 
latter was obtained from the @-component of the 
equation of motion. The pressure distribution is 
obtained from the corresponding r-component of 
the equation of motion : 


dp /dr (25) 


p (v,?/r) 


This equation could be integrated directly to 
obtain the radial pressure distribution at a given 
elevation if the pressure were given at some value 
of the radial co-ordinate r. In viscometry, 
however, situations may arise in which one 


knows only the pressure py of the fluid at rest at 
the given elevation, and that the volume and 
fluid content of the annulus are unaffected by 
rotation. The latter problem is the one considered 
here. 

The problem may be solved by temporarily 
relaxing the assumption of incompressibility. It 
can be shown from equation (2) and (3) in com- 
ponent form that equation (4) with the left-hand 
side equal to zero, and equation (25), are still 
valid for the isothermal steady flow of a compres- 
sible fluid provided that one assumes that » does 
not depend upon the pressure. These equations 
then have to be solved along with the equation 
of state which is here taken to be the first two 
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Fic. 1(a). Steady-state pressure profiles when both Fic. 1 (b). Steady-state pressure profiles when the inner 
cylinders are moving with the same angular velocity. cylinder is stationary. 
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Fic. 1(c). Steady-state pressure profiles when the outer 
cylinder is stationary. 
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terms of a Taylor series for p as a function of the 
Gibbs free energy G = H — TS: 


Po + ~ Go) 


(sa), (¢ 


po [1 + 6 (G — G,)] 


p 


mii +d | pap) 
r. 
in which p, and p, are the uniform pressure and 
density of the fluid when the system is at rest, 
and b, is the value of b= (1 p)(Q9 IG)y— (Vp W)y 
at p Pe 
The 


performed by inserting the velocity distribution 


integration of equation (25) may be 
of equation (10) and assuming that p is known as 


a function of p: 


P 

1 2 

p ‘dp v|r(s 
i 


P 


in which N and @ are 


dimensionless quantities : 


es 


*)@)" 
1 @ 


(R AQP AB 


(28) 


B 
A 


The constant of integration C is determined by 


the normalizing condition that 


Po £ dé pé dé 


which ts just a statement of conservation of mass. 
For the equation of state in equation (26) this 
reduces to: 


1 p 


| | p dp Ed€ 0 


< Pe 


Substitution of equation (27) into equation (31) 


vields the integration constant C and hence the 
final expression for the pressure distribution : 


P 


“ey 


3(é @2 z2iné 


p ! dp N | 


Ps 
= 
—- In «| (32) 


«x 


For incompressible fluids the left side of equation 


Po), Po- 
has been obtained by expressing the equation of 


(32) just becomes (p The same result 
state as a Taylor expansion of p in powers of 
(p — po). but the development is not as simple 
as that given above. Three important limiting 


cases for the pressure distribution are ; 
(i) Q: 


Pp ~ Po 


l 
pal? R 7," + «?) (33) 


Q, = 0, 


2, # 0 (or « =0 and Q = x): 
P—Po\(' — **) 

pylde? 4 | - 

1, 
(' 


« 


2iné 


8 x? 


In - 
+ x? 1 ne? 


0 (ors 


land@ — 1): 


Zin é 


+ 2 x? 
In x 
an 


(35) 


These functions are shown graphically in Fig. 1; 
clearly the largest pressure gradients correspond 
to equation (33). These results may be used to 
determine whether or not pressure effects are a 
disturbing factor in any given concentric cylinder 
viscometric measurements, They may also be of 
use in calculating radial pressure distributions in 
lubrication systems. 
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NOTATION 


functions defined in equation (10) 

quantity deferred in equation (26) 

static pressure 

radial co-ordinate 

radius of outer cylinder 

tangential velocity 

dimensionless angular velocity (equation 8) 
ratio of radius of inner cylinder to that of outer cylinder 
coefficient of shear viscosity 

dimensionless radial co-ordinate (equation 5) 
fluid density 


dimensionless tangential velocity (equation 7) 
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Precise testing of binary vapour-liquid equilibrium data by the 
Gibbs=—Duhem equation 


H. C. Van Ness* 


Department of Chemical Engineering, Rensselaer Polytechnic Institute, Troy, New York 
(Received 12 March 1959) 


The general methods of testing thermodynamic data for binary solutions by the 
test is 


Abstract 
Gibbs-Duhem equation are briefly reviewed, and a new “* composition—resolution ~ 


proposed. Thermodynamically exact methods of testing binary vapour-liquid equilibrium 
data, both for constant pressure and for constant temperature, are then developed. Application 
of the 


* composition-resolution ” test is fully discussed, and its advantages are demonstrated. 


Résumé 


solutions binaires par léquation de Gibbs-Duhem ont été résumées et une vérification nouvelle 


Les méthodes générales de verification des données thermodynamiques pour des 


* composition-résolution “ a été proposée. Les auteurs ont développé des méthodes thermody- 
namiquement exactes du contréle des données d’équilibre d'un binaire vapeur liquide a pression 
constante et aussi a température constante, L’application de la vérification composition 
résolution est complétement discutée et ses avantages sont démontrés, 


Zusammenfassung Die allgemeinen Methoden zur Priifung thermodynamischer Daten fiir 
biniire Lésungen durch die Gibbs Duhem Gleichung werden kurz mitgeteilt und eine neue 
Ferner werden thermodynamisch exakte Methoden zur 


Zerlegungsmethode vorgeschlagen. 


Priifung der biniren Gleichgewichtsdaten dampfformigtlissig, und zwar sowohl fiir konstanten 


Druck wie fir konstante Temperatur entwickelt 


Die Anwendung der Zerlegungsmethode wird 


ausfibrlich mitgeteilt und ihre Vortecile aufgezeigt. 


Data for solutions are most conveniently dealt 
with by putting them in the form of excess 
property changes of mixing. If G is taken to 
the excess 


represent a general, 


properties of mixing are given by : 


AG 


property in 


AG* Ages’ 


AG® — AG 


A (ite! 


where AG“ represents the property change 
which would result if an 
formed, and the bar indicates a partial molal 
property change of mixing, i.e. AG; — G; — G,. 

There are three general methods for testing the 


ideal solution were 


thermodynamic consistency of experimentally 
determined partial molal excess properties of 
binary systems. The first is known as the * 
test,” and has been recognized for many years. 
For data at constant temperature and pressure, 
it requires that : 


* slope 


P © 464) . = AG5) 


(3) 


da, das 


In other words, for thermodynamically consistent 
data a definite relationship must exist between 
the slopes of the curves of AG" and AGE vs. a ,. 
The disadvantage of this test is that slopes are 
very difficult to determine accurately, and the 
uncertainty of the method is usually at least as 
great as the uncertainty of the data. 

The second method is known as the “ area 
test.” A simple and general derivation is given 
by Van Ness and Mrazex [1]. For data at 
constant temperature and pressure it requires 


that : 
1 1 


| AGE dx, = [ AGg ar, (4) 


That is, the areas under the curves of AG® and 


* During 1958-59, Visiting Fulbright Lecturer in Chemical Engineering at King’s College, University of Durham, 


Newcastle upon Tyne. 
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AG§ vs. x, should be equal. This is an integral 
test, and suffers from the disadvantage that it 
does not represent a point-by-point test of the 
data. It is entirely possible for this integral or 
“area” test to be satisfied when the “ slope 
test ” at various concentrations is not. 

The third method might be called the ** compo- 
test. It has apparently not 
been discussed before, and represents a new 
general method for the testing of data by the 
Gibbs-Duhem equation which is straightforward 
and accurate. Furthermore, it allows a point-by- 
point evaluation of the data. It is based on the 
observation that a given curve for the excess 
property change of composition 
(v4 Or @,) can result from an infinite number of 
sets of experimentally determined values of AG* 
and AGE, but only one set of values of AG* and 
AG* can result from a given curve for the excess 
In other words, if 


sition-resolution *’ 


mixing vs. 


property change of mixing. 
the experimental values of AG and AG) are 
combined according to the equation : 


AG® — a2, AG" +2, AGE (5) 


and the resulting curve of AG* vs. 2, is then used 
(for example, by the method of tangent intercepts) 
to calculate values of AG" and AG4,, the derived 
set of values should agree with the experimentally 
measured set. If they do, the data are consistent. 
If the data are of poor quality, this method can 
do no more than reject them. But if they are 
reasonably accurate, the result of this process of 
testing will be a set of data smoothed in such a 
way as to be completely consistent. 

The “composition” step of this process as 
represented by equation (5) presents no difliculty,. 
But the “ resolution ” 
classical method of tangent intercepts is usually 
not sufficiently precise for effective application 
of the method. Van Ness and Mrazex [1] have 
developed a new method of tangent intercepts 
which results in the accurate determination of 
partial molal excess properties. It is based on a 
plot of AG® /r, a, vs. #4, or alternatively on a 
reciprocal plot of a, 2,/AG* vs. v4. Many of 
the excess property changes of mixing which 
AG*® may represent are set forth in the paper 
cited [1]. Of these, the one most conveniently 


step carried out by the 


calculated from vapour—liquid equilibrum data 
is AF® ‘RT, where AF® is the excess Gibbs free 
energy of mixing of the liquid phase. Its relation 
to the activity coefficients of the liquid phase is 
given by: 

AF® /RT 


vain ¥4 +2gln Vp (6) 


Phase equilibrium data for binary systems 
cannot be taken with both temperature and 
pressure held constant, as has been assumed in 
the foregoing discussion, and the method of testing 
is not quite so direct. Consideration has to be 
given to two cases: constant-temperature data 
and constant-pressure data. In either case the 
general relationship given by equation (6) is 
valid. The other exact equation which relates 
the activity coefticients and AF®/RT for any 
phase has been derived by the author [2] : 


AH 


d1 ? ~ y, Ad; 7 
RT? iT (In y;da,) (7) 


‘E : 
d | AF At dp 


RT RT 
where AF*®, Av, AH, and the ¥;'s are all taken 


with respect to the pure components at the same 
temperature the 
physical state as the phase considered. We will 
apply this equation to a binary liquid phase in 
equilibrium with its vapour, first, at constant 


and pressure and in same 


temperature, and second, at constant pressure. 
For a binary liquid at constant temperature, 

equation (7) reduces to : 

d(AF® RT) 


dx, 


In ¥, —In ¥p 4 (=) (3) (3) 


Total derivatives may be used for the case con- 


sidered because there is but one independent 
variable, here taken as 2,. Fig. 1 shows a typical 
plot of AF*®/RT vs. x, for constant pressure. 
The slope of a tangent drawn to the curve is, of 
course, d(AF® RT)/dx,, and the tangent is 
shown intersecting the ordinates at 7, = 1 and 
2,4 = 0 at points labelled « and £ respectively. 
If the data for AF* RT were for constant 
temperature and pressure, then « and 8 would 
be partial molal properties and indeed would be 
identical with In ¥4 and In Y,. But in the case 
considered where the pressure varies, « and £8 
differ from In Y, and In ¥, by amounts which are 
determined as follows: « and £8 are related to 
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Constant 7 or constont 0 











"4 


Fic. 1. 
constant-temperature or constant-pressure data. 


Classical method of tangent intercepts applied to 


AF*® RT by equations which are obvious from 
inspection of Fig. 1: 
AF® = 
RT . B 
AF® 


P= ee ~ “a 


d(AF*® RT) 
da A 
d(AF*/RT) 


da i 


(9a) 


(9b) 


If equations (9a) and (9b) are each combined 
with equations (6) and (8), the resulting express- 


7) (a) 


d 
7) (a,) 


and 


ions for « and 8 are: 


In Y4 (10a) 


+ a, (= 
2(2 


exact, 


In ¥p (10b) 


show the 


These 


difference between the tangent intercepts « and 8 


equations are 


and In Y, and In Y, for a plot of AF* RT vs. a4 


at constant temperature. In_ practice, this 
difference is negligible, because Av is 
almost always very small. Hence it 


entirely satisfactory to treat the tangent intercepts 


usually 
is usually 


as being identical to the In ¥;'s. 

For the case of constant-pressure data, equation 
(7) becomes : 
d(AF*®/RT) 

dx 4 


A plot of AF*, RT vs. x, for data at constant 


= Invy,—In7,- 


AH \ , dT 
R rT) | 


) (11) 


dx, 


Dg % 


VaN NEss 


pressure is identical in form to the plot of constant- 
temperature data, and is also represented by 
Fig. 1. Indeed, equations (6) and (9) are equally 
valid for this case, 
equations (10a) and (10b) are readily derived 
by combining equations (9a) and (9b) with 
equations (6) and (11); 


and equations analogous to 


InY¥, —tp 


aaa) (5) (12a) 


RT? 
AH \ (dT 
Rr) (i..) 


Again we note the distinction 
intercepts, « and 8, and Iny, and In ¥,. 


B=In Yp, (12b) 
between the 
Here 
this difference cannot usually be neglected, and 
it is not in general satisfactory to regard the 
intercepts as representing values of InY, and 
In? »- 

In either case, constant temperature or constant 
the treatment of data is in principle 
Equations (10) and (12) can both be 


pressure, 
the same. 
written : 


In v4 (13a) 


In yg (13b) 


l 
(=. (Z,.) for constant-temperature 
data, 


“— (omy (at 


and 6 RT? et constant-pressure 


dx 4 


It should be noted that : 


+ tpB 


(However, « 4 In ¥, and 8 # In ¥,). 


“EE 


RT 


rginY¥, +apln7p, (14) 


The intercepts, 
equation : 


x and £, are related by the 


*) 


which follows from equations (9a) and (9b) after 
differentiation of « with respect to 2, and of 8 
with respect to vy. Equation (15) indicates the 
appropriate of phase-equilibrium 
data calculated from experimental measurements 
by equations (13a) and (13b). These same data 


(i, 


dx, 


“slope test ” 
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may be subjected to an 
to the equation : 


‘area test’ according 


1 1 


[ «ar, | Bda, 


. 
v0 0 


(a — B)dr1, = 0 


0 


(16b) 


But a- B In V4 = In YB + 0 


Therefore : 
Y 
In“ da, 


‘B 
0 0 


dda, (17) 


This particular test was discussed by THiJssen 
[3]. Both the “ 
suffer from already 
The “ composition—resolution ” test proposed here 


slope test ” and the “ area test ” 
disadvantages discussed. 
is also applicable, and represents a much more 
satisfactory procedure. 

Values of AF*® RT are calculated from the 
experimental data, and are plotted vs. 2 ,. Values 
of « and 8 can then be determined by the usual 
method of tangent intercepts, and the results are 
used to calculate values of In Y, and Inv, by 
equations (13a) and (13b). If the derived values 
of Iny, and In ¥, correlate with the values 
calculated directly from the data, then the data 
are consistent. 

As previously explained, there is considerable 
inherent uncertainty in drawing tangents to the 
curves of AF®/RT vs. 2, and in extrapolating 
to get values of « and 8. The procedure for 
overcoming this problem has been described by 
Van Ness and Mrazex [1]. Rather than plotting 
AF® RT vs. x, one plots AF* x,a, RT vs. x4 
as shown in Fig. 2. The tangent intercepts on 
this plot, 7, and /,, are related to « and 8 by the 
equations ; 


(18a) 


AFE , 
“a4apy RT ° 


, APE 
2 ry RT | 


B - (18b) 


and In ¥4 and In Y, are given by : 
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Duhem equation 
AF® 
42%, RT 


InYz v P AFE 
Ly, 41" 242, RT 


Iny 
“4 = rp |2 
TR 





Constant 7 or constont p 








Pe se ee ee ee 





A4 


Fic. 2. Method of Van Ness and Mrazek applied to 
constant-temperature or constant-pressure data. 


It should be noted that: 


AF* 
v42, RT 


Inv, , n¥p 


(20) 
Lp v4 
Thus the comparison here is most conveniently 
made between experimental and derived values 
of (In ¥4) 2, and (In ¥,)/ 7,4. 
superiority of this procedure are given in the 


The reasons for the 


paper cited [1]. It is also explained that in some 
cases the reciprocal plot of 7,2, RT /AF® vs. x, 
can be used to greater advantage. 

For constant-pressure data the rigorous testing 
procedure requires a knowledge of the heat of 
mixing as a function of composition and tempera- 
ture. In practice this means the determination 
of the heat of mixing as a function of composition 
at two temperatures. Since the heat of mixing 
at a given composition is almost always essentially 
linear with temperature, interpolation for the 
effect of temperature is very easy. Unfortunately, 
at present the necessary heat-of-mixing data are 
usually lacking. 

For constant-temperature data it is usually 
satisfactory to take 6 as zero. Hence the recom- 
mended testing procedure requires no additional 
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data beyond phase-equilibrium measurements. 
In Fig. 3 data of Réck and Scurdper [4] for 
the acetone-chloroform system at the constant 
temperature of 30°C are subjected to the * compo- 
sition-resolution " test. The 
corrected for the non-ideality of the vapour 
This appears to be essential even for 


data have been 
phase. 


sub-atmospheric pressures, because the method 





Zz On IF /x A, RT 











Fig. 3. Acetone — chloroform system at 30 © [4). 


is very sensitive to deviations of any kind. The 
points shown on Fig. 3 represent values calculated 
from the experimental data. The 
AF* (x42, RT vs. x4 has been drawn through 
the data points so as to give a smooth relation. 
The curves of (In ¥,)/a, and (In ¥4)/24 vs. @4 have 
been derived from the AF* 2,2, RT curve by 
equations (19a) and (19b), with 6 taken as being 
negligible. These curves are seen to correlate the 
experimental points reasonably well, and it is 
concluded that the data are essentially consistent, 


curve of 


Moreover, the derived curves represent the data 
smoothed in such a way as to be completely 
consistent. Curves of the smoothed values of 
Iny, and In¥, are also shown on Fig. 3. It 
should be noted that if the curve of AF* ‘x, 2, RT 
is expressed in equation form, analytical express- 
ions for Iny, and In¥, can be readily derived 
from it, If this curve is a straight line, then the 
2-constant Margules equations result. If the 
reciprocal plot, 2,7, RT/AF® vs. x4, is a 
straight line, then the Van Laar equations are 
obtained. 

For constant-pressure data, where 5 cannot 
usually be neglected, the quantities derived from 
the AF® xv, 2, RT vs. x, curve are x and 8B (see 
equations 18a and I8b] and not Iny, and 
In ¥,-. Hence derived analytical expressions are 
for « and 8, and again Margules and van Laar 
equations arise (for « and §) if the appropriate 
curves are straight lines. Correct equations for 
In ¥, and In ¥, can be obtained only by expressing 
5 as well in analytical form. In general, the 
nature of 5 is such that this is most difficult, and 
it is probably best treated graphically. 

An alternative procedure for the case of 
constant-pressure data is to make use of AF*® 
itself rather than AF* RT. A plot of AF* vs. 2, 
at constant pressure would be similar to Fig. 1. 
In this case we will designate the tangent inter- 
cepts by «’ and 8’. The required equations are : 


d AF® ) 


21 
dx 4 ( * 


AF*® 4 ry ( 


(21b) 


if ae 
p’ = AF* o.(—— ) 


dx, 


AF*® vr, RT WY, +2,p RT iny, (22) 
By differentiation : 


d(AF®/RT) 


dx 4 


or (s-) 1 (" = 


RT? \dx,) | RT\ dv, 


If this is combined with equation (11), the result 
is: 

i AF*®  (/AF* — AH) (dT 
oof (05 iS). 

dix, T di, 
+ RT \Iny, — RT ny, 
or 
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d AF® 


dx A 


RT ny, —RTiny, — AS*® (=~) (23) 
A 


Combination of equations (21), (22), and (23) 
gives : 

dT 

». 
dz) (24a) 


=) 


dx 4 


RT In¥4 —*y As* ( 


B= RT ny, +24 AS® | (24b) 
In practice, of course, one should plot AF” /x, x, 
vs. v4 and use the method of Van Ness and 
Mrazek for determining values of x’ and 8’. 

A possible advantage of equations (24) over 
equations (12) is that AS® may be small even 
though AH is not. This is the basis for the many 
suggestions in the literature that it is preferable 
to deal with T In y,; rather than with In 7, itself 
for constant-pressure data. Although this is not 
necessarily true, it is sometimes the case. It is 
doubtful, however, that the AS* term of equation 
(24) is very often negligible. 

Another possible advantage to the use of AS” 
rather than AH is that AS* may not vary nearly 
so rapidly with temperature. In this case values 
of AS* need be determined at just one tempera- 
ture, whereas AH must be known as a function 
of temperature. 


The area test in terms of these variables is 


readily deduced from either equation (23) or 
(24), and requires : 


| RT in "4 dx, | ase (“7 ) ar, 
Yr da, 
0 0 

Thus it is seen that the suggestion of BrouGuTon 
and Brear.ey [5] that the left-hand integral in 
equation (25) is zero is not in general correct. 

One other question often arises in connexion 
with constant-pressure data: Do the curves of 
In Y¥, and In Y, vs. 2, have to be tangent to the 
horizontal at their respective origins ’ The 
answer is an unequivocal yes. This may be seen 
by differentiation of equation (13a), which gives : 


din ¥, d 


dir, 


dx, 


As 2,4 approaches unity, zp, 5, and da/dz, all 


Duhem equation 


approach zero, and all terms on the right-hand 
side of this equation That 
approaches zero can be seen from equation (15). 


vanish. da dx, 
This conclusion is based on the premise that none 
of the derivatives, dx dx,, d 8 dx,, or dT dx ,, 
becomes infinite as 2, approaches zero or unity. 

At an azeotropic point 6 becomes zero for both 
constant-pressure and constant-temperature data. 
Moreover, the curves of AF*® RT vs. a, for 
constant pressure and for constant temperature 
which pass through a given azeotropic point are 
tangent at this point. 

It has been the object of this paper to develop 
thermodynamically exact equations for the treat- 
ment of vapour-liquid equilibrium data at both 
constant pressure and constant temperature. In 
addition, precise methods for the testing and 
smoothing of such data have been presented. 


NOTATION 
excess Gibbs free energy change of mixing 
a thermodynamic property in general 
a molal property of pure i at the T and p of 
the solution 
a partial molal property 
a property change of mixing per mole of 
mixture 
a partial molal property change of mixing 
an excess property change of mixing 
an excess partial molal property change fo 
mixing 
a property change of mixing for an ideal 
solution 
a partial molal property change of mixing for 
an ideal solution 
heat of mixing per mole of mixture = Al” 
tangent intercepts as shown in Fig. 2 
pressure 
universal gas constant 
excess entropy change of mixing per mole of 
mixture 
absolute temperature 
volume change of mixing per mole of mixture 
mole fraction of component i in solution 
tangent intercepts at 2 4 = 1 as shown in Fig. 1 
tangent intercepts at 24 = 0 as shown in Fig. 1 
- deviation term defined following equation (13) 


activity coefficient of component i in solution 
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Abstract 


for the case when Dorcy’s law is valid and the Péclet number is large. 


This paper presents the results of an analysis of dispersion in a porous medium, 


It is shown that, after a 


sufficiently long time, the effective longitudinal diffusivity is a logarithmic function of the Péclet 


number. 


Résume 


Dorey s‘applique et pour des nombres de Péclet considérables. 


Analyse de la dispersion d'un écoulement dans un milieu poreux, lorsque la loi de 


On montre qu’apres un temps 


suffisant, la diffusivité longitudinale est une fonction logarithmique du nombre de Péclet. 


Zusammenfassung 


Untersuchungsergebnisse 


tiber die Durchmischung in einem pordsen 


Medium werden mitgeteilt im Bereich des Doreyschen Gesetzes und bei grosser Pécletzahl. Nach 
geniigend langer Zeit ist die effektive Langsdurchmischung eine logarithmische Funktion der 


Pécletzahl. 


1. Description or THE Mopet 
One. of the difficulties in calculating the properties 
of flow through porous media is that the detailed 
geometry of the medium is not usually known 
and the flow field cannot be computed in detail. 
Theoretical work is at present confined to the 
investigation of models which can be handled 
mathematically. This note summarizes an investi- 
gation (SarrMan [2]) of the dispersion of a 
dynamically neutral material quantity in flow 
through a particular model, first considered by 
JosseLin pe Jone [1], and suggested independ- 
ently by Sir Grorrrey TayLor. 

The model consists of a random, statistically 
isotropic, network of straight circular capillaries 
of length | and radius a, several capillaries starting 
and finishing at each junction. The case is con- 
sidered in which the Reynolds number of the 
flow through the capillaries is small so that 
Poiseuille flow exists in each capillary or pore. 
To a good approximation, the pressure gradient 
along a pore making an angle @ with the direction 
of the mean flow can be taken as cos @ dP dW, 
where ) P Xr is the mean pressure gradient 
and the mean flow is supposed uniform and in 
the a-direction. The velocity profile in this 
capillary is 


l MP 
. 2 2 
q by: cos @ = (a r?) (1) 
where r is the distance from the axis of the pore 
and y« ts the viscosity of the fluid. The proportion 
of pores with @ lying between @ and @ + d@ is 
sin 0d @, and the mean interstitial velocity (in 
the x-direction) is therefore 
jn a 

¢ 2r 7 )P , 
| | q cos 0 ; drsn@dé@ A l 

4 a 


. Say (2) 
24u Wr : 


0 0 
The filter velocity is o U', where o is the porosity, 
and therefore the permeability of the model is 


I 
24 o a’ (3) 


In terms of U’, (1) may be written 


q = 6U cos@(1 r? a*) (4) 


LONGITUDINAL DISPERSION 
Fiutp PARTICLE 


Tue 

OF A 

Consider the motion of a fluid particle. After 
passage through one capillary it has gone a 
distance | cos@ in the -direction in a_ time 
t = lq. Since the network is random, it may be 
assumed that the values of @ and r for successive 
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pores traversed by the particle are random 
variables and _ statistically independent. The 
probability distribution of the displacement and 
time after a fluid particle has passed through n 
pores can then be calculated, using well known 
methods of statistics. It is then possible to obtain 
the probability distribution (p.d.) of the longitu- 
dinal displacement, 2, of the particle after a 
given time T. 

It is found that the mean displacement is given 


r UT (5) 


by 


so that the centre of a cloud of marked fluid 
particles moves with the interstitial velocity U. 

For the dispersion about the mean, it is found 
that (2 — ¥)% 
= — o, and that 


is asymptotically normal as 


(a 7)? 547 (6) 


l 


l 
~ lo 
tS | 8 


Ir 


where the bar denotes a mean, and it is necessary 
that log (%1) >» 1 before equation (6) applies. 

It is to be noted that the p.d. of 2 — # is not 
asymptotically normal ; it has a tail for negative 
values of a 
of the probability distribution is negligible when 
weighted with 7" ?, but not otherwise. 

If an effective diffusivity, k,, for the longi- 
tudinal dispersion is defined by the formula 
(x — ¥)* = 2k, T, then k, Ul is not constant but 
is an increasing function of ¥ /. This result is not 


%. whose effect on the cumulants 


in accordance with diffusion processes in general, 
and arises here because the time of passage of a 
fluid particle through a pore, ¢t, becomes infinite 
as r +a or 0 (1 2). If the quantity being 
dispersed is subject to molecular diffusion with 
a diffusivity D, it is clearly not physically sensible 
to suppose that elements of the material quantity 
may spend an infinite time in a pore, however 
small D may be. It is therefore necessary to 
investigate the effects of molecular diffusion, even 
though D <« k,. 


3. Errects OF MOLECULAR 
DIFFUSION 
Molecular diffusion acts in two ways: (i) it 
transports material across a capillary, the time 
scale for this process being a*/8D = t,, say (cf. 


TayLor [3]) ; and (ii) it transports material along 
a capillary with a time scale [/2D = t,, say. 
These processes are negligible if ¢, and t, are very 
much greater than lg, where q is to be interpreted 
as the velocity of the fluid particle at the point 
where a material element enters a capillary. But 
however small D, provided it is not identically 
zero, there are values of r and @ for which these 
conditions are not satisfied. 

An exact calculation of the effects of molecular 
diffusion presents great difliculties. However, if 


>t, > UU (7) 


molecular diffusion can be taken into account in 
an approximate manner by assuming that its 
effect is to impose “ cut-offs’ on the possible 
of t. Two “cut-offs” are required: 
firstly to express the idea that material near the 
capillary walls will diffuse away from the walls 
in a time of order ¢, to where it will be convected 
with a velocity of order 3U cos @, and secondly 
that material diffuses along a pore in a time of 
These may be expressed formally by 
saying that the time spent by a material element 
in a capillary is the least of l/q, t, + 1/3U cos @, 
and ¢,. One reason for believing that this rough 
approximation is not too unreasonable lies in the 
fact that the results depend logarithmically upon 
D and are not sensitive to variations of t, and t,. 

The p.d. of the displacement of material 
elements after a given time can then be calculated 
and it is found that the mean displacement still 
moves with the velocity U, i.e. equation (5) is 
satisfied. The dispersion about the mean now 
depends, however, on the relative values of the 
dimensionless quantities #1, Ut,/l and Ut,/1. 
The following results are obtained, where the 
square of the dispersion is denoted by s? = (« — ¥)*, 
which are valid when (7) is satisfied. 

(a) If, roughly, (%/1)"? > Ut, l > Ut,/1 > 1, then 


- (log sa ) 


values 


order t,. 


og Uo 
l 


a more precise condition for the validity of (8a) 
being 
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Ut,/l 


(#1)? (log 3Ut,/1)"* * (8b) 


7/1 = N, say, is proportional to the average 
number of pores which a fluid particle goes 
through in the time 7. 1/t, = 2D/l =, say, 
can be called the * velocity of diffusion along a 
pore,” and //t, = 8Dl/a* = v,, say, the “ velocity 
of diffusion across the pore.” Then Ut, 1 = Ul 2D 

U vland Ut, 1 = Ua? /8DIl = U v,. The con- 
dition (7) can be expressed as U >» v, > v,, and 
the rough condition for case (a) as 
N¥?2 > Uv, > U/v,. This can be satisfied by 
taking N sufficiently large, however small D may 
be. 

Expressions for other relative values of the 
parameters are : 


If, 


‘ 


Ut/l > (¥/))"? > Ut,/l > 1, 


(1 


(b) roughly, 


1 


12 


more precise conditions being 
7 (7) > , 3Ut, /l 
BN) > @1)*flog(z 1/2}? * 

3Ut,/l 


(#/1)* *{log(z/1)¥2})2 (9b) 


> 1 
(c) If, roughly, Ut,/l > Ut,/l > (%, 1)? > 1, then 

2 
(log ; 


more precise conditions being 


7\'" Ut, 1 
log | i) rs G8 log & 
4Ut,/l 


(z/l)’? log (¥ 


s? 


zl 


1 
. 48 
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j (10a) 


[jh 2 ” 


pila > (10b) 


It may be that ¢, and ¢, differ by a factor of 
100 or more, in which case the condition 


(11) 


may possibly be satisfied, i.e. the time scale for 
diffusion of material across the capillary is small 


t> l//U>t, 


-/ 


compared with the time spent in the pore. In 
this case, the material will be uniformly concen- 
trated across a cross-section and the material is 
effectively convected along the pore with the 
average velocity 3U’ cos@ (see Taytor [3)). 
The analysis leading to the above results can then 
be repeated taking the time spent by a material 
element in a capillary as /1/(3U cos @), and 
assuming that the effect of molecular diffusion 
is to “ cut-off ” the possible values of the time by 
an upper bound ¢,. The results valid when (11) 
is satisfied are that ¥ = UT and as follows: 


(d) _' > 1, then 


If, roughly, (%/1)? » =a 


gs? 
zl 


= *s 3Ut, 1 
at ee 12 


more precise conditions being 


Ut,/l 


(z/l)2 (log 3 "t, (12b) 


Jyre2 <1 


(e) If, roughly, Ut,/l > (%/1)"? > 1, then 


1 277 
oO ss 
6 e 2l 


s2 


= * 
zl aan 


more precise conditions being 


3Ut, 1 


(F 1)? flog (1)! 2}1/2 > 1. (18b) 


log (7) > 1 


The precise magnitude of in these 
formulae is not certain, but interpreting it as 
“ greater than by a factor of five” should not 
lead to serious error. 

The p.d. of (a — %)/%"? is in each case approxi- 
mately normal, but that of « — % is not, having 
a long tail for negative values of « —¥. The 
third moment has been calculated for cases (a) 
and (d), and it is found that 


(@ =a ~ — Se Ult, > coast + wo. (14) 


(The third moment of (# — ¥) 7? is O(#-?) 
which + 0 as ¥ — o). 


* This is effectively the result given by JossELIN DE 
Jone [1)}. 
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4. Discussion 

Although the theory applies strictly only to a 
network of capillaries, it is expected that it will 
have some relevance, qualitatively at least, to 
dispersion in flow through an actual porous 
medium, with the length of the capillaries, /, 
identified with the grain size of the medium, and 
the radius a estimated from the permeability 
and porosity by means of equation (3). It is to 
that the 
quantitatively accurate because they depend upon 


be noticed, anyway, results are not 
the precise form of the cut-off conditions, although 
not to any great extent. 

For given values of U and D, and ¥ varying, 
the dispersion must eventually be given by cases 
(a) or (d) since the conditions for one of these 
cases can always be satisfied by taking 7 suflfici- 
ently large ; however, the larger U'] D, the larger 
must be # The dispersion is then proportional 
to # 2, in accordance with the usual behaviour of 
diffusion logarithmic 


function of Ul D. (Compare Taytor [3] where 


processes, and is also a 
dispersion in a straight capillary is found to be 
inversely proportional to D). 

If D 0 (ie. t and U 


kept fixed, case (c) eventually obtains, which is 


and t, + «), with # 


a 
of course identical with (6) which was derived 
(10b) therefore 
molecular 


on the assumption that D — 0. 
the that 
should be negligible. Notice that log ( /)' * must 


gives condition diffusion 
be large compared with unity for the distribution 
to be gaussian, which implies values of ¥ / of at 
least 2 « 104 
proportional to 7? log ¥, although in practice the 
difficult to 


In this case the dispersion is 


logarithmic dependence would be 
observe. 

Equivalently, suppose that 7 and PD are fixed 
and U varies. For the small (but not too small) 
values of U, cases (a) or (d) will obtain. As U 
increases, the dispersion increases, because, in 
effect, molecular diffusion has less time to decrease 
concentration gradients, and (c) eventually 
obtains. 

These comments are illustrated in the figure, 
where s* %/ is plotted against Ul D on a logarith- 
10° 


The full curves are the values 


mic seale, with, for the sake of example, % / 
1/5. 
given by the various results in the region of their 


and al 




















109. UO 


Fig. 1. 
s* rl as a function of Ul D, for #/1 = 108 and al 
The full lines give the results for the case of the accompany- 
ing letter, and the dotted lines are obtained by interpolation. 
The horizontal dotted line is the value given by case (e). 


validity and the letter denotes the corresponding 
case. It will be noticed that much of the range 
is not covered by the cases, but the value of 
s? 7l can then be obtained by interpolation, 


shown by the dotted lines on the figure. The 


values shown for very large values of Ul/D may 


be academic, since the value of U’ cannot be too 
large else Poiseuille flow will not exist in the 
Roughly this requires 
1, where v is the kinematic viscosity of 


capillaries. speaking, 
lay 
the thuid. 
given by (13a), Le. the dispersion according to 
(e) if the conditions for that case were satisfied. 
One point seems worth emphasizing, and that 
is that the results of laboratory experiments with 
relatively small values of % 1 must apparently 
be applied with care to situations occurring in 
practice where the value of ¥ 1 may be very 
different, especially if Ul D is extremely large. 


The horizontal dotted line is the value 


5. Tur LATERAL 

It is found that the lateral or transverse 
displacement in a direction perpendicular to the 
mean motion has a p.d. which is normally distri- 
buted with zero mean and variance. 


y—3 8UlT 


provided (%/1)"? » 1. 


DISPERSION 


3,8 lz, (15) 


This expression is equiva- 
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lent to an effective lateral diffusivity k, — (3. 16)1 2. 
This result is independent of molecular diffusion, 
provided of course D << Ul. However, its deriva- 
tion assumes that the direction of successive pores 
through which a marked fluid particle moves are 
statistically independent, and there are reasons 
based on continuity arguments for believing that 
this assumption is not valid for the lateral 
dispersion, although reasonable for the longitu- 


dinal dispersion. The lateral dispersion may 


therefore be an order of magnitude less than that 
given by (15). 

A full discussion of the assumptions and a 
description of the analysis are given in [2]. Also, 
the theory is there compared with some experi- 
mental data, and the agreement is found to be 
not unreasonable. 


NOTATION 
radius of capillary 


molecular diffusivity 
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effective longitudinal diffusivity 
effective lateral diffusivity 
permeability 

length of capillary 

mean pressure 

velocity in a capillary 

distance from axis of a capillary 
longitudinal dispersion 

cut-off time for diffusion across a 
capillary 

cut-off time for diffusion along a 
eapillary 

axial component of the interstitial 
velocity 

mean longitudinal displacement in 
time T. 

viscosity of fluid 

kinematic- viscosity 

porosity 

angle between directions of capillary 
and the mean flow. 
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Abstract 


Mass-transfer rates in a disk-column were studied by vaporizing liquids in a range 


of counter-current flowing gases. The results of these studies are well correlated by a “ relative 
jp factor,” which takes into account the relative velocity of the gas and liquid streams. 
The importance of the gas-liquid relative velocity in the analysis of packed-column absorp- 
tion studies is illustrated by showing how the true effect of diffusivity can be masked. 
For a disk-column, a simple correlation is established between the gas-film mass-transfer 
coefficient, the diffusion coeflicient and the gas-phase pressure drop. The relationship provides 
excellent support for the 2/3 exponent on the diffusivity. 


Résumé 


liquides se vaporisant dans un gaz circulant a contre-courant. 
reliés convenablement par un “ facteur relatif jp’ 


courants de gaz et de liquide. 


L’auteur étudie les vitesses d'échange massique dans une colonne a disque avec des 


La résultats de ces études sont 


qui tient compte de la vitesse relative des 


L’importance de la vitesse relative gaz liquide, de 'absorption dans les colonnes de garnissage 
est mise en évidence en montrant combien leffet exact de la diffusion peut étre masque. 

L’auteur établit pour une colonne a disque une relation simple entre le coefficient de transfert 
de masse gauz-film, le coefficient de diffusion et la perte de charge dans la phase gazeuse. La 
relation montre la valeur de lexposant 2.3 du coefficient de diffusion. 





fas ung 


Zusamn 


von Fliissigkeiten in entgegenstromende Gase untersucht. 


In einer Scheibenkolonne wurden Stoffiibergangswerte durch Verdampfung 


Die Ergebnisse lassen sich gut durch 


einen “ relativen jp-Faktor  wiedergeben, der die relative Geschwindigkeit zwischen Gas und 


Flissigkeit beriicksichtigt. 


Die Bedeutung dieser relativen Geschwindigkeit bei der Untersuchung von Fillkérper- 
Absorptionskolonnen wird dadurch bewiesen, das gezeigt wird, wie die wahre Wirkung der 


Diffusion iiberdeckt werden kann. 


Fiir eine Scheibenkolonne wird cine einfache Beziehung zwischen dem Stoffiibergangskoe- 
flizient auf der Gasseite, dem Diffusionskoeflizient und dem gasseitigen Druckabfall angegeben. 
Der Exponent 2 3 fiir den Diffusionsterm wird gut bestitigt. 


INTRODUCTION 
Tue disk-column, as introduced by Sreruens 
and Mornis [1], has been used primarily for the 
investigation of gas absorption systems in which 
the liquid-film is the major diffusional resistance. 
All previous measurements of gas-film controlled 
mass-transfer rates in a standard disk-column 
{1}, [2] and [3] have been made with air as the 
principal gas phase. Further investigation was 
considered desirable to establish the effects of the 
properties of other gases on the gas-phase mass- 
transfer process. It seemed reasonable to expect 
that such a study might also elucidate some of the 


factors generally influencing mass transfer: e.g. 


the relative velocity of the fluid streams; the 
true effects of diffusivity and the relation between 
the pressure drop and the mass-transfer coefficient. 
It is possible to investigate these variables in a 
disk-column, because of the known gas-liquid 
interfacial area. 

In this investigation, the gas-phase mass- 
transfer process was studied by measuring the 
rates of vaporization of liquids in different gases 
for a wide range of counter-current gas and liquid 
flow rates. Using this technique, complications 
arising from diffusional resistance in the liquid 
phase were absent and the gas-film mass-transfer 
coefficients could be measured directly. 
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EXPERIMENTAL 
During a preliminary investigation, the liquid 
velocity on the disks was determined by hold-up 
For 
water, at liquid flow rates less than 240 lb/hr ft 
these velocities compared favourably with those 


measurements on an eighty disk assembly. 


calculated according to the method of Coorrr, 
et al. [4]. At greater liquid flow rates, a departure 
from the theoretical equation was evident. The 
observed velocities were less than those obtained 
by calculation. This is presumably due to the 
formation of ripples on the rims of the disks at 
high liquid rates. Ripple formation was reported 
by Taytor and Ropers [3]. 

In the experiments, the gas-phase mass-transfer 
process was studied by measuring the rates of 
vaporization of water into counter-current streams 
of hydrogen, air, carbon dioxide, dichlorodifluoro- 
methane and the vaporization of n-heptane by 
air in a standard disk-column. The column and 
the auxilliary equipment used has been described 
previously by the author [5]. The principal 
dimensions of the column are given in Table 1. 
The column was operated under conditions such 
that no spray or splashing of the liquid was 
observable. To achieve these conditions in the 
n-heptane experiments, extremely small liquid 


rates were required. 
Table 1. dimensions 


Disk-column, principal 


Number of disks | 27 

O-161 ft? 
O-O4S4 ft 
O-O148 ft 
0-os2 ft 
89°, 
0-122ft 
O-O52 ft 





Dry surface area of disks* 
Average disk diameter 

Average disk thickness 

Tube internal diameter 

Free space (dry) 

Mean perimeter for liquid flow 
Equivalent diameter for gas flow 





*Uncorrected for liquid film thickness or loss of area at 
points of contact, but used throughout as an approxima- 
tion of the true liquid surface area. 


The partial pressure of the diffusing component 
in the outlet gas stream was measured con- 
tinuously with a calibrated gas analysis cell of 
the thermal conductivity type. The inlet gas to 
the column was analysed at the beginning and 


end of each series of run. The change was too 


The 
equilibrium vapour pressures were interpolated 
from the extensive tables compiled by Strut [6). 


small to warrant continuous analysis. 


The diffusion coefficients for the water vapour 
systems were taken from the table of selected 
experimental data given by WILKE and Lee [7] 
and the diffusion coefficient for the air-n-heptane 
system was obtained from the measurements of 
SCHLINGER et al. [8]. The relevant diffusion data 
are given in Table 2. 

The rates of vaporization were calculated from 
the steady-state partial pressure of vapour in the 
exit gas stream, the inlet partial pressure and the 
gas flow rate. For the n-heptane experiments, 
the calculated true mean driving force was used 
the 
However, in the water vaporizations, because of 


to evaluate mass-transfer coeflicient. 
the very much higher liquid rates, negligible 
error was introduced by using the log mean of 
the terminal driving forces, based on bulk liquid 
The calculated 
transfer coetlicients and other relevant results for 


temperatures. gas-film mass- 


the systems studied are given in Table 3. 


Table 2. Published diffusion data 





Diffusion 
coefficient 


(em? sec) 


Schmidt 
number 


System Temp. 
(C) 





0-850 
0-260 
O-164 
0-105 
0-073 


H,-H,O 
Air H,O 
CO,-H,0 
CCLF,-H,0 
Air n Cry, 


kwew nw 
or Ge Ge 





CORRELATION DISCUSSION 
or RESULTS 


AND 


For a disk-column, SterHens and Morris [1] 
and Taytor and Roperts [3] have shown that 
the relative velocity of the gas and liquid streams 
can be used to correlate results of countercurrent 
and co-current absorptions of ammonia in water 
from a dilute ammonia-air mixture. The basis 
for using the relative velocity of the gas and 
liquid streams for a particular type of absorption, 
i.e. either countercurrent or co-current has not 
been investigated before. 
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Table 3. Experimental and calculated results 





Partial pressure 
Liquid | Gas Liquid Relative of diffusing Temperature Gas film 
rate velocity | velocity velocity component (°C) | mass transfer 
(Ib, hr ft) (ft sec) (ft sec) (ft sec) (mm Hg) ' am -| coefficient 
— + —| Liquid | Liquid Gas Gas Ib mole 


- - 
Inlet | Outlet inlet | outlet inlet outlet (hr ft? atm) 


‘ — me: ‘a 

Vaporization of water with hydrogen 

2: O+ | 138 | 25.4 24.4 26.3 
O53 | 20° 23-8 26-2 
0-5 244 | 235 | 26-2 
0-6 25- 23-1 25-9 
OF 22-8 25-4 
O-4 23-6 25-9 
O-4 230 | 260 
Os 228 | 26-2 
OS 220 26-0 
O-6 21-5 25-0 
0-6 220 23-8 
7 22:! 23-6 


oO? 23- 


= 
= 


1.23 
1-70 
2-038 
2-26 
2-42 
1-35 
2-05 
2-26 
3-02 
3-30 
2-76 
2-39 
1-86 
1-29 
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te 


te 
oe 


os 
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.- .- - - 
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> Ww 
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te nw 


tS ts 


7 24-8 


Vaporization of water with air 
1 24-0 1 

2% | 240 a 

12° 24-0 23-6 23° 1: 

1. 


2 
eo 
— =F ot St St of Ge UF 


= 


te t te 


13: 24-0 14 
14 240 om 
14 23-9 35 2: O75 
16- 23-0 2: 2 O-67 
18-6 25-8 O-60 


' te t te 
te 
= 
=~ 
- 


26-2 
26-1 
25-8 


25-3 


0-60 
0-59 
O58 


O52 
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: a a a 5 
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O74 
OTs 
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0-62 
O97 
O-s9 
O-S4 
O-s80 
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Table 3—contd. Experimental and calculated results. 





Partial pressure 
of diffusing 


component 


Relative 
velocity 
(ft sec) 


Liquid Gas 
rate | 


(Ub he ft) | 


Liquid Gas film 
| mass transfer 
coefficient 
Ib mole 
(hr ft? atm) 


Temperature 
(C) 
(mm Hg) - t — 
| -_— ' Liquid 
Inlet Outlet inlet 


velocity 
(ft (sec) 


velocity 

(ft /sec) 
Gas 

outlet 


Gas 


Liquid 
outlet inlet 
Vaporization of water with carbon dioxide 

6 | 188 26-7 26-2 
18-4 26-9 26-3 
16-7 | 269 | 263 
14-2 26-9 260 | 
16-4 aL 26-2 
3-55 17-1 270) 26-4 
317 17-5 7 26-5 
2-73 18-4 7 26-6 

a) 


2-59 29-4 
2)1 
29-0 


26 


200 
3-14 
75 


30 


-—-— = 


27: 


to to 


7 
5 
28-3 
29-0 
20-9 
31-0 
31° 


~ 


te t 


2-16 3 19-8 | 26- 


1-06 21- 71 26-4 


Vaporization of water with dichlorodifluoromethane 
1-69 0 15-8 23:1 | 22- 
zol 0 13-7 23-1 
2-22 0 13-4 23-1 
15-6 23-1 
12:5 23-1 
15-1 23-1 
14-4 23-1 


~ 


~ 
oo 


=. 
- 
~ 


tw Ff 


wo Se 


219 
264 
1-92 


1-87 0 


os os 
~ © 


R 
a) 
5 
w25 
. 
7 
5 


te te t 


Vaporization of n-heptane with air 
1-83 0 8-0 13-7 f 12: 
2-65 0 6-1 13-1 35 | 32s 
285 | Oo 48 12: 12: 
540 | O 4-2 12° 11-3 
3:8 12-2 11-7 


6-90 0 
RRS 0 3-2 12-4 11-7 


15 
15 
15 
15 
15 


2 
15 


~ 
=> 
te 


Ss = 

. a 

te te 
ow ot Se oe oe oe 


0-25¢ 
0-30t 
0-39 
O-45F 
0-48¢ 


-~ 
= 
te 


- 
= 
te 


-~ 
= 
te 





*These temperatures measured close to the last disk by copper-constantan thermocouple with cold junction at 0°C, 
+Corrected for true mean driving force. 


In an earlier publication [5] the author has 
shown that a conventional Reynolds number vs. 
friction factor type of graph is not satisfactory 
for correlating the gas-phase pressure drops 
across a disk-column, when the disks are irrigated 
with water. A further analysis of these results 
has since been made [9], in which the significant 
velocity term was taken as the relative velocity 
of the gas an liquid streams. Pressure drops across 
the dry column and those when irrigated with 
water at rates less than that required for ripple 
formation can all be correlated by a single line 


on a “relative friction factor” vs. “ relative 


Reynolds number” graph. These findings were 
considered to support the contention that the 
relative velocity of the gas and liquid streams is 
also the significant velocity term in the mass- 
transfer process. 

If the gas-liquid relative velocity were signi- 
ficant in the original dimensional analysis of the 
mass-transfer process, then a modified approach 
could be derived. 

(1) 


? (=?) (5): ( 


A “relative j,” factor can then be defined : 


kL 


1.e. - D 
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kK. Pas { # x (2 


Relative j 
- UR P oD 


{ Ip) 


. 23 
k, af He i (3) 
Trp Pe ob 

It is proposed that the “ relative j,, factor ~ 
should be a function of the “ relative Reynolds 


number 


(Rip 


(fn). (4) 


ya 


The newly detined “ relative j, factor” and 


the conventional j,, correlation were applied to 
the experimental results, Fig. 1 is a comparison 


of the two methods of correlation, The relative 


Cia “Tio FT 


Reictive Reynolds number 


Comparison of the j,, and (jp)py correlations, 
(a) Constant liquid rate I $25 lb he ft 
Water vaporization 

hydrogen 
@ ar 
D carbon dioxide 
dichloroditluoromet hane 
oO air [ 208 lb brft 
(b) All liquid rates. 
Water vaporization 
hydrogen 
@® «air 
BD carbon dioxide 
dichlorodifluoromet hane 
n-Heptane vaporization 
air 


The 


jp” correlation includes all liquid rates. 


jp correlation is only for the constant liquid rate 


of 325 Ib hr ft. 

The conventional j,, type of correlation shows 
the hydrogen results as being displaced below 
the air points, particularly at low Reynolds 
numbers, where the effects of relative velocity 
are greatest. The difference between the air and 
hydrogen poimts decreases with increased Re. 
This is to be expected, if the relative velocity 
is significant, as at higher values of Re, the liquid 
velocity has a much smaller effect, because of 
The dichloroditluoro- 
the 


results 


the higher pas velocities, 


methane points, however, do not exhibit 


same trend. To be consistent, these 
should be displaced above the air whereas they 
are also displaced below. 

rate, the newly 
in vs. “relative Re” 


plot correlates all the results except the dichloro- 


Irrespective of liquuid flow 
proposed =“ relative 
difluoromethane points, which are displaced well 
below the line of best fit through the other points. 
A similar discrepancy is observed for the two 
These two 
bulk of the 


represent cases 


low velocity carbon-dioxide points. 
the 
points 


carbon-dioxide points and 
dichlorodifluoromethane 
in which the liquid velocity on the disks is greater 
than the upward gas velocity. 

At low gas velocities, the gas is carried down- 
ward by the liquid and the resulting reverse 
flow mixes the gas vertically in the column so that 
than for true 


the mean driving force is less 


countercurrent contact. Under these conditions, 
mass-transfer coetlicients, evaluated by assuming 
no longitudinal mixing, are considerably smaller 
the Coorer et al. [10] 


reported similar anomalies in a packed tower 


than true coetlicients. 


when the downward liquid velocity exceeded the 
upward gas velocity. 


VELOCITY 
COLUMNS 


RELATIVE 
PACKED 


LiguID 
EFFECTS IN 


Gas 


The effect of diffusivity on the rate of gas-phase 
mass-transfer in packed-columns has been the 
subject of 


contlicting — reports 


anomalies for many years. 


and apparent 
In general, it would 
appear the experimental results for systems in 
which there is no liquid flow (e.g. pre-wetted 
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raschig rings, naphthalene berl saddles, ete.) 
can be well correlated using the j, factor with 
the 2,3 the Schmidt number, 


whereas for irrigated systems a range of exponents 


exponent on 


from 0-15 to 0-90 has been required to correlate 
the results. Suutman et al. (11) has shown that 
differences in effective interfacial areas can 
explain the low exponent of 0-15 obtained by 
Surosky and Dopce [12]. The exponents on 
the Scumipr number of 0-90 and 0-77 obtained 
by Lyxcn and Wintke [13] and Yosnuipa [14], 
respectively, cannot be explained by differences 
areas, one liquid 


in interfacial because only 


phase was used. It is proposed that the true 
effect of diffusivity in these studies is masked 
by the relative velocity of the gas and liquid 
streams. 

Although a formal analysis of relative velocity 
effects in a packed-column requires knowledge 
of the flow paths of both the gas and liquid, 
the effect of the gas liquid relative velocity on the 
exponent of the Schmidt number can be demon- 
strated by evaluating the relative velocities in 
terms of an unknown tortuosity (8) of the gas 
flow path and an assumed true mean liquid 
velocity (v,). The true mean gas velocity can be 


calculated in terms of 8 by equation (5). 


3(“). 


pe 


Now since Ug 
(6) 


For example, the results of Lyxcu and WILKE 
[13] can be correlated by the “ relative j, factor” 
with the 2.3 exponent on the Schmidt number, 
if the magnitude of v, 8 in equation (6) is assumed 
to be 0-1 ft see. This assumption corresponds 
with a liquid velocity of 0-067 ft ‘sec when the 
value of 8 is 1:5. Both of these figures are of an 
order which would be expected in a_packed- 
column, The results are 
plotted as the product 8 [a (jp)g] vs. 1/8 (vq p/#) 
on a log-log scale in Fig. 2. The points plotted 
cover the range of Lyxcu and WILKe’s experi- 
ments up to the observed loading points. 

The good correlation of the points in Fig. 2 


selected recalculated 











Fig. 2. Correlation of Lyxcu and Witke [13] data using 
‘relative jp factor” and * 


Data of Lyxncu and Wiik1 


relative Reynolds number.” 
liquid rate 1575 Ib hr ft® 


O-1 ft see. water vaporization 


t 
recalculated assuming —£ 
~ 


ws 
A helium 
@ air 
dichlorodifttuoromet hane 


shows very clearly the marked effect of the small 
liquid velocity correction. In the main, this is 
due to the very low apparent dichlorodifluoro- 
methane 0-25 to 1-4 ft/sec. 


The exponent of the Schmidt number in this 


gas velocities of 
correlation is 2.3 as compared with the exponent 
of 0-9 obtained by Lyncu and WILKE [13] when 
using the conventional correlation. 


Mass-TRANSFER 
THE 


rHeE 
TERMS OF 


Drop 


The gas-phase pressure drop across the disk- 


CORRELATION OF 
CORFFICIENTS IN 

PRESSURE 
column was measured with a micromanometer 
for each of the mass transfer runs. These readings 
were converted to the pressure drop across one 
disk by correcting for the pressure drop due to 
the column shell and by dividing by the number 
of disks in the assembly. A formal theoretical 
analogy between the pressure drops and the mass- 
transfer coetlicients was not possible because of 
the unknown contribution of form drag to the 
total drop. <An_ entirely 


empirical approach was used to determine if any 


measured pressure 
consistent relationship existed. 

For a given system, it seemed reasonable to 
expect the pressure drop and the mass transfer 
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coefficient to be affected similarly by changes in 
gas velocity and liquid rate. A plot of k, versus 
AP, for the air-water system showed that these 
quantities could be correlated by a straight line 
on a log-log graph as shown in Fig. 3. This method 
of correlation is very attractive, particularly 
as the uncertain effects of liquid flowing over 
the disks are eliminated, because the mass-transfer 
coeflicient and the pressure drop are affected 
equally by these uncertainties. 





Gas tiim coefficient 





5 4 
SP, 1 'C 


3. The relation between pressure drop and the 


mass-transfer coefficient. Air 
Water vaporization with air 
I’ — 825 lb/hr ft 


water system. 


208 Ib her ft 
240 lb hr ft 
140 Ib her ft 
64 Ib he ft 





Gas film coetticen 10° 





4P, 1 \0* 


Fic. 4. Generalized correlation of gas-film mass-transfer 
coeflicients in terms of the diffusivity and the pressure drop. 
Water vaporization 
A hydrogen 
@ air 
© carbon dioxide 
+ dichlorodifluoromethane 
n—Heptane vaporization 


@ air 


All the mass transfer results of this investiga- 
tion can be well correlated by a log-log plot of 
k, D-** vs, AP, as shown in Fig. 4. This type of 
correlation has not been reported previously. 
The 2/3 exponent on the diffusivity is considered 
to lend excellent support for the use of this 
exponent in correlating gas-phase mass-transfer 
results, because the uncertain effects of liquid 
velocity, associated with other published investi- 
gations, are eliminated in this approach. 
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NOTATION 
area of interphase contact in a packed 
ft? fr 
ft? hr 
Ib hr ft® 


column 
molecular diffusivity 
gas superficial mass velocity 
ke Pam { # )" . 
Ug Pm oD 


dimensionless 


hy PRM | # ) . dimensionless 
UR Pm eb 
mass-transfer coeflicient 

Ib mole, hr ft? Ib mole /ft® 
mass-transfer coeflicient Ib mole /hr ft? atm 
characteristic dimension ft 
atm 
pressure drop across a disk Ib ft? 
logarithmic mean of the partial pressure of 
the phase 


total pressure 


the non-diffusing component at 


boundary and the bulk of the fluid atm 


Lv p 
Reynolds number — dimensionless 


“ 
Lup p 
“ 
dimensionless 
ft /hr 


= “relative Reynolds number ™ 


fluid velocity 
average gas velocity ft /hr 
average liquid velocity ft /hr 
relative velocity of the gas and liquid streams 
ft/hr 
tortuosity of the gas fluid path in a 
dimensionless 
dimensionless 


packed column 
packed bed porosity 
I’ = liquid rate per wetted perimeter 
Ib /hr ft 
Ib /ft® 
Ib moles /ft® 
Ib /ft hr 


p = gas density 
Pm = molar gas density 
ps = gas viscosity 
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Chemical Engineering Practice. 
Cremer and T. Davies. Vol. V. 
Butterworths Scientific Publications, 
Academic Press, New York, 1958. vi 
$13.30 


Systems I. 
and 
XXiV pp. 
London 95s. 


Chemical Engineering Practice. Hdited by H. W. 
Cremer and T. Davirs. Vol. VI. Fluid Systems II. 
Butterworths Scientific Publications, London and Academic 


Press, New York, 1958. 600 pp. 95s. $13.30. 


Ir is fashionable at the present time to regard the * unit 
operations * concept in chemical engineering as obsoles- 
cent, and this idea is certainly being fostered in the 
chemical engineering schools. Nevertheless, there are no 
chemical engineers who do not know what is meant by 
unit operations, and it is with certain of these transfer 
processes, to use the current description, that most of 
5 and 6 of Chemical Engineering Practice are 
books titles Fluid 
Systems I and II respectively. 

The first half of Volume 5 follows logically from Volume 


4, which contained a large section on the academic ideas 


Volumes 


concerned, although the carry the 


of fluid mechanies. In Volume 5, these ideas are supple- 
mented by a useful practical discussion of the selection 
and performance of pumps, fans, blowers and compressors, 
and descriptions of pipe systems for liquid and gases. 
Worked examples of design calculations are included 
for those who regard this series as handbooks of procedures. 
Next in the book is an authoritative chapter on the design 
of high pressure vessels, in which the academic origin is 
very obvious, but which does contain a lot of useful 
practical information. The next chapter on the production 
of vacuum combines a good description of the various 
kinds of equipment with a valuable section containing 
worked examples on flow rates and pumping rates for the 
molecular flow regions, pump-down times and so on. 
The rest of Volume 5 consists of a 327 page contribution 
All the 


usual text-book material is included in a commendably 


on solvent extraction and distillation processes. 


brief form, supplemented by critical comments on the 
latest work The 
chapters in this section have the headings, “ Liquid 
liquid equilibrium ; 
extraction equipment ; vapour liquid equilibria ; 
tion of binary mixtures ; 
and distillation equipment.” The authors are all associated 
with the Department of Chemical Engineering at Birming- 
ham and are known for their original contributions to 


and by comprehensive references, 


liquid-liquid = stagewise contact ; 
rectifica- 


multi-component distillation, 


the subject. The practical equipment side is adequately 
covered in the space available. Taking the section as a 
whole, it appears that this will be the starting point of a 
literature survey for many future students. 


Volume 6 continues the discussion of gas-liquid systems 
by chapters on gas absorption and on the operating 


research 


columns. Once 
again, there is the apparently mixture of 


text-book and a good deal of practical information and 


characteristics of packed absorption 


inevitable 


design data obtained from actual operating plant, although 
The 


final chapter dealing with gas-liquid systems is a very 


much of the latter has been previously published. 


good one on evaporation practice, with the emphasis 
The of text and 
illustrations is very informative indeed. 

next. 


really on the practice. combination 


Gas- solid considered There are 


chapters on fluidization and fluidized beds, and on the 


systems are 


general applications of fluidization with special reference 
to the carrying out of chemical reactions. These sections 
contain an orthodox discussion of fluidizing velocity, 
pressure drop, mass and heat transfer and a description 
of some of the processes for which fluidized beds are used, 
together with some good practical points. Anyone ex- 
pecting up-to-the-minute information will be disappointed 
at the omission of such topics as the effect of bubble 
recently attracted a deal 


The other topic which can be considered 


formation which has good 
of attention. 
under the heading of gas- solid systems is Adsorption and 
the chapter on this topic is a general description of the 
elementary principles, together with useful illustrations 
of how the process is used in industry. 

Chapters on the liquefaction and fractionation of gases 
and on refrigeration practice can best be described as a 
useful introduction to these subjects, and are accompanied 

The next main 
is mainly concerned with the 


by exhaustive references. heading, 
* Liquid-liquid Systems ~ 
separation processes of leaching, crystallization, filtration 
and centrifuging. The authors of these chapters have 
contrived to blend an adequate amount of the simple 
theory with a gocd deal of useful practical information, 
particularly in the selection and operation of filters. There 
is a chapter on colloids, which comprises nearly all a stu- 
dent needs to know in a degree course, and the Volume 
is completed with a chapter on sublimation and vacuum 
freeze drying. 

If an attempt is to be made to assess the value of the 
work as a whole, reference must be made to what 
H. W. Cremer, the General Editor, wrote in the preface 
to Volume I of the series, “ The chemical engineer is 
increasingly at a loss to find ready-made answers to his 
problems. There still remains a great gulf fixed between 
theory and practice, a space which, of necessity, he must 
on occasion himself endeavour to fill.” Do these volumes 
help him? The blend of practical information with the 
academic approach which is characteristic in all the con- 
tents certainly deserves to be successful. And perhaps 
Mr. Cremer would that since he wrote his 
preface the gulf between theory and practice is narrowing 


perceptibly. 


now agree 


E. S. SELLERS 
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Low Temperature Physics and Chemistry. Edited 
by J. R. Dituincer. University of Wisconsin Press, 
1958. 676 pp. $6.00. 


Tuts book contains summaries of the 225 papers presented 
at the Fifth International Conference on Low Temperature 
held at the University of Wisconsin in August, 1957. 
It is obviously impossible to review such a book but it 
may be useful to the reader to list the main headings 
under which the papers were grouped. These follow : 
I. Viscosity and flow properties of liquid He?. 
Il. Heat flow in, and specific heat of, liquid He?. 
III. Dielectric constant. Refractive index, 
thickness and rotation of liquid He*. 
IV. Heat diffraction, 
structure of He, 
Liquid He*® and He. 
Theory of liquid He* and liquid He®. 
Solutions of He® in He*. 


sound transmission and 


The temperature scale and temperature measure- 
ment 
heat 


Instrumentation: He*® Dewar 


liquetier, 


systems, 

switches, a computer elements ; 
techniques. 

Theories of superconductivity 

Properties of some superconductors. 

Intermediate state of superconductors. 

Superconductivity of thin films, contact resis- 
tance, pressure effects and volume changes in 
superconductivity. 

Millimeter and the 


magnetic effect in superconductors. 


wave absorption para- 
Thermal and electrical conductivity of elements 
and compounds. 
Thermal and electrical conductivity of alloys. 
Measurements of specific heats of solids. 
Analyses of specific heat data. 
De Haas-Van Alphen and other oscillatory 
effects. 
Magnetoresistance, thermoelectricity, resistivity 
and hall effect. 
Mechanical properties of solids including 
solidified gases 
Paramagnetism and resonances. 
Magnetic susceptibilities. 
Ferromagnetism and general magnetic properties. 
Thermal expansion, absorption spectra, dielectric 
properties, rotational transitions, polymers 
and nuclear magnetic resonance. 
XXVI. 
It can be seen from this list that the papers deal almost 
entirely with phenomena at the temperature of liquid 
helium and below. There is practically nothing which is 
of direct interest to the chemical engineer. 
B. F. Dopcr 


Nuclear orientation experiments. 


O. Kupascuewski and E. Lil. Evans: Metallurgical 
Thermochemistry. (Third Edition — revised). Pergamon 
Press, London, 1958, xiv + 426 pp. 63s. 


Reavers of the earlier editions, and in the field of physical 
chemistry there must be many indeed, will recall that 
the book was divided into five sections, namely: The 
The Estima- 
Thermochemical Data ; 


Theoretical Basis ; Experimental Methods ; 
tion of Thermochemical Data ; 
Examples of Thermochemical Treatment of Metallurgical 
These third 


The first and third sections remain materially 


Problems. divisions are retained in the 
edition. 
the same in this later edition, whilst in the second section 
there have been a few changes illustrating new apparatus 
and techniques, particularly in the part dealing with 
calorimetry. ‘Two additional examples in the fifth section 
are to be welcomed, and the tables of thermochemical 
data in the fourth section have been revised and extended. 
It is disappointing, however, in this very useful fourth 
section, to find that the data on dilute solutions, which 
included infermation on the solubilities of gases in metals, 
have been omitted from this edition. 

One cannot help feeling that anyone already possessing 
a copy of the second edition would have perhaps preferred 
a supplementary volume containing only the revised and 
extended tables of thermochemical data, particularly 
as the third edition of the book has increased in price by 
eight shillings. 

However, for those interested in modern thermo- 
chemistry, this book remains one of the outstanding works 
in its field. One can argue the theoretical section would 
have to be supplemented by additional reading, and that 
the experimental section is uneven in that the review of 
calorimetric methods is more complete than, for example, 
the parts dealing with cells and gas-condensed phase 
equilibria measurements ; it is certain, though, that the 
last three sections would immeasurably benefit anyone 
who has to handle thermochemical data. 

The book obviously has a metallurgical bias but the 
writer believes that this in no detracts from its 
general appeal. It is the emphasis on the critical assess- 
ment of thermochemical data and the correct manipulation 
of such data, together with the reference tables, which 
make the book so valuable to anyone who is interested 
in the practical application of chemical thermodynamics. 


way 


B. C. H. STEELE 


Electrochemical Processes in Chemical 
1957. 464 pp. 30s. 


\. REGNER: 
Industries. Constable. 


On account of the title one would suppose that all elec- 
finding application in industry 
This however is not the case: the book 


trochemical processes 


are covered. 


describes only processes in aqueous solutions and so the 
electrolysis of molten salts (manufacture of aluminium, 
sodium, ete.) and the electrochemistry of gases (e.g. the 
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manufacture of ozone) are not dealt with. On the other 
hand some processes, such as the manufacture of hydro- 
chloric acid from salt and sulphuric acid, which have no 
relation at all to electrochemistry, are included. 

The first 170 pages contain a survey of theoretical 
electrochemistry. The objection could be made that this 
part is rather superfluous as, in the first place, there are 


already many books covering this subject and, secondly, 


because many of the theoretical explanations are not 
referred to when dealing with the applications. So, for 
instance, the theory of the conductance of strong elec- 
trolytes, that of galvanic cells and that of the cathodic 
precipitation of metals which have no follow-up in the 
second part. 

The description of many processes is more elaborate 
than may be found in other books covering the same 
field (Binuarer, MANTELL) which naturally is an advantage. 
On the other hand however obsolete processes are treated, 
e.g. the old processes for the manufacture of bleaching 
powder, the electrolytical manufacture of hypochlorite 
and the Griesheim-Elektron cell which are of no more 
than historical interest. 

The bibliography is not at all adequate. All that is 
given is a list of about thirty books, only twelve of them 
being less than twenty years old, and half of these are in 
Russian. One of the consequences is that one always 
doubts whether the given process data originate from 
own experience or are taken from literature. The same 
applies for some explanations of reaction mechanisms 
which are sometimes not very likely. For instance, is it 
not true that the anodic evolution of chlorine from a 
chloride solution is due to overvoltage of oxygen, neither 
is the reasoning why ammonium sulphate gives a higher 
output of persulphate than sulphuric acid as an electrolyte 
very convincing. 

In conclusion it must be said that this book has hardly 
any advantage above the existing textbooks and owing 
to its lack of literature references, it is useless as a reference 
book. 

J. G. 


HooGLAND 


W. MiaLk1: Kernverfahrenstechnik. Eine Einfihrung 
fiir Ingenieure. Springer-Verlag, Berlin/Géttingen/ 
Heidelberg 1958. 4835S. 78,60 DM. 


Es gibt in der Bundesrepublik noch nicht viele gute 
Einfiihrungen in die Reaktortechnik. So hat der Verfasser 
eine echte Liicke ausgefiillt, indem er in Anlehnung an die 
vorziiglichen amerikanischen Einfiihrungen VON 
GLASSTONE, STEPHENSON und Murray ein deutsches 
Lehrbuch verfasste. Aus der speziellen Fachrichtung des 
Verfassers ist es wohl zu erkliren, dass er sein Buch 
* Kern-Verfahrenstechnik * genannt hat, obwohl es sich 
mehr um “ Reaktortechnik " handelt. Beim Bemihen, 
die Verfahrenstechnik dennoch nicht zu kurz kommen zu 
lassen, mussten im Stoff Kompromisse geschlossen werden, 
die der Einheitlichkeit und Ausgewogenheit des Buches 
abtriglich sind. 

Dass der Verdampfungsvorgang bewusst nicht beriick- 
sichtigt wurde, ist angesichts der Wichtigkeit der Siede- 
wasserreaktoren allzu einseitig ; demgegeniiber wurde der 
Wiarmeiibergang bei fliissigen Metallen eher zu ausfiihrlich 
behandelt. Die Moderatoren sollten wegen ihrer Wichtig- 
keit eingehender besprochen, die Kiihlmittel nicht in der 
Gliederung als Werkstoffe bezeichnet werden. 

In einer Neuauflage sollte cine Reihe von Fehlern 
(z.B. Coulomb’sches Gesetz, 5. 83-87) oder Ungenauigk- 
eiten (z.B. Diffusionslinge L und Fermialter 7 sind als 
wirkliche mittlere Entfernungen angegeben) vermieden 
werden. Auch sollten mdglichst fir die Wirkungsquer- 
schnitte die allgemein tiblichen Bezeichnungen gewahit 
werden, also 


“total 
7 Absorption 


"Absorption T %Streuung ; 
*Spaltung * °Strahlungseinfang - 


Zu den unbestreithbaren Vorziigen des Buches gehdrt 
das wohlgegliickte Bemiihen des Verfassers, die physi- 
kalischen Vorgiinge anschaulich zu behandeln, wodurch 
es sich fiir Studierende empfiehlt. 

K. JAroscHek 
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Errata 
T. J. Gunsert: Liquid mixing on bubble-cap and sieve plates. Chem. Engng. Sci, 1959 10 243. 


Equation (4) on p. 249 should read as follows : 
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N. L. Karencavz and T. D. Bocnaroyv: Effect of height of a flat rectangular conduit on heat transfer to water. 
Teploenergetika 1959 6 (3) 76-78. 





